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ABSTRACT 


In an effort to reduce weight while maintaining high strength, many contemporary 
structural systems are designed with lower margins of safety than predecessors. The 
criterion for minimum weight design is predominant in the design of aircraft, missile, and 
spacecraft vehicles. One obvious means of obtaining a high strength, minimum weight 
design is the use of light, thin-walled structural members of high strength alloys. Thin- 
walled beams of open sections such as I, Z, Channel and angle sections are frequently 
used for intricate structures in spacecrafts. Due to low torsional rigidity of thin-walled 
beams of open sections, the problem of torsional vibrations and stability is of prime 
interest. For the last three to four decades mechanical vibrations have been recognized as 
a major factor in the design. Mechanical vibrations produce increased stress, energy loss 
and noise that should be considered in the design stages if these undesirable effects are to 
be avoided, or minimised. 

The research work carried out particularly deals with dynamic analysis of lengthy 
uniform thin-walled composite I-beams of open section. Fibre-Reinforced Plastics (FRP) 
have been increasingly used over the past few decades in a variety of structures that 
require high ratio of stiffness and strength to weight. The analytical model developed 
by Lee and Kim [11] is taken and the dynamic behaviour of a thin-walled I-section 
composite beam is studied in detail. This model accounts for the coupling of flexural and 
torsional modes for arbitrary laminate stacking sequence configuration, i.e. symmetric as 
well as unsymmetric, and various boundary conditions are also discussed in detail. A 
displacement-based one-dimensional finite element model is developed to predict natural 
frequencies and corresponding vibration modes for a thin-walled composite beam and the 
equations of motion are derived from Hamilton's principle. 



Further experiments are conducted by taking two specimens of glass-epoxy 
composite rectangular laminated beams of different fibre orientations and natural 
frequencies are determined for a cantilever beam by impact experiment. These results are 
compared with the natural frequencies of composite rectangular laminated beams by 
using ANSYS R10.0 .the numerical results from finite element analysis showed in 
general a good agreement with experimental values. 

The present research work is extended with the assumption that the finite element 
analysis of laminated composite I-beams will be valid and the effect of fibre angle, 
modulus ratio, height-to-thickness ratio, and boundary conditions on the vibration 
frequencies and mode shapes of the composite I-beams is studied. Using ANSYS R10.0, 
general purpose FEA software) the modal and harmonic analysis are carried out and the 
results obtained are presented in graphical form, for uniform and tapered composite I- 
beams. The modal and harmonic analysis is also carried out for a cantilever composite I- 
beam for two cases, with its bigger end fixed and smaller end free and again with its 
smaller end fixed and bigger end free. The study was extended to know the behaviour of 
the same composite I-beam with circular holes in the web. 
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Chapter 1 


INTRODUCTION 

1.1. Introduction to thin- walled beams 

For intricate structures such as space craft, beams of standard cross sections may not 
be the most efficient or convenient structural members to use. Thin walled beams of open 
sections are frequently employed for their structural efficiency. When the dimensions of cross 
sections of a beam are very small compared to length of the beam, the beam is referred as 
Thin-walled beam. Thin- walled beams of open sections such as I, Z, Channel and angle 
sections are frequently used for intricate structures in spacecrafts. The thin-walled beams of 
open sections are widely used in engineering industry and space vehicles due to their low 
torsional rigidity and stability. With the improved versions of extrusion methods in metal 
forming, beams of different shapes of cross sections can be formed to order. However, it has 
been known that a beam with nonsymmetrical cross sections under loads will not only deflect 
but also twist. Only under special loading along the flexure axis, beam will deflect without 
twist. 


Currently, fibre-reinforced laminated composite structures are being used extensively 
in a variety of industrial products and structures such as air craft, automobiles, marine 
vessels, rocket launchers, pressure vessels and pipes owing to their strength, high stiffness, 
low weight and good resistance to fatigue and corrosion. The most common problem is 
delamination along ply interfaces which may arise due to fabrication defects and impact 
loading. The present work particularly deals with dynamic analysis of lengthy uniform thin- 
walled I-beams of open section, particularly composite beams. In the present thesis, the 
analytical model developed by Lee and Kim [11] is taken as reference and the dynamic 
behaviour of a thin-walled I-section composite beam is studied in detail. This model accounts 
for the coupling of flexural and torsional modes for arbitrary laminate stacking sequence 
configuration, i.e. symmetric as well as unsymmetric, and various boundary conditions are 
also discussed in detail. 

Some of the following basic assumptions made in this study are 

• The thin-walled beam has uniform open cross section along its length. 



• Strains are assumed to remain within the elastic limit and the curvature and twist 
of the beam are considered to be small. 

• The material of the beam is homogeneous and isotropic and obeys Hook’s Law as in 
the case of conventional beams. 

• The centroid and shear centre of the cross-section coincide. 

• Shearing strains of the middle surface due to shear and warping effects and axial 
strains are considered to be negligibly small. 

• Longitudinal inertia effects are considered to be very small. 

• Damping forces are not considered. 

• Distortion of cross section in their own planes is not considered. 

The need for high performance to weight ratio structures in the most advanced 
engineering applications leads to massive usage of composite materials for critical 
applications. The thin-walled composite I- beams are modelled and modal and harmonic 
analysis is done to know the effects of fibre angle, modulus ratio, height-to-thickness ratio 
on the vibration frequencies and mode shapes of the composites by using ANSYS, general 
purpose software. The different terminology used in analysis software is discussed below. 

1.2. Convergence Criterion 

The solution should converge to the correct result, if the following three conditions are 
achieved or satisfied. 

• The displacement function must be continuous within the element choosing 
polynomial for displacement model. 

• The displacement function must be capable of representing rigid body displacement 
of the element i.e. When the nodes are given such displacement corresponds to a rigid 
body motion, the element should not experience any strain and hence leads to zero 
nodal forces. The constant terms in polynomials used for the displacement models 
would usually ensure this condition. 

• The displacement function must be capable of representing constant strain states 
within the element when the body or structure is divided to smaller and smaller 
element. As the number of elements increases to infinity the strains in each element 
also approach to constant value. 



Hence the assumed displacement function should include terms for representing constant 
strain states. For ID, 2D, and 3D elasticity problems, the linear terms present in the 
polynomial satisfy the requirement. 

1.3 Criterion of Finite Element Selection 

Element to be used is dependent on physical problem itself. Following are important 
considerations in the selection of type of element. 

1. Number of degrees of freedom required. 

2. Expected accuracy. 

3. Easy in formulating the necessary equation. 

4. Degree to which the physical structures can be modelled without approximation. 

1.4 Defining the Layered Configuration 

The most important characteristic of a composite material is its layered configuration. 
Each layer may be made of a different orthotropic material and may have its principal 
directions oriented differently. For laminated composites, the fibre directions determine layer 
orientation. 

Two methods are available to define the layered configuration: 

1. By specifying individual layer properties 

2. By defining constitutive matrices that relate generalized forces and moments to 
generalized strains and curvatures (available only for SOLID46 and SHELL99 ) 



1.4.1 Specifying Individual Layer Properties 


With this method, the layer configuration is defined layer-by-layer from bottom to 
top. The bottom layer is designated as layer 1, and additional layers are stacked from bottom 
to top in the positive Z (normal) direction of the element coordinate system. 

Occasionally, a physical layer will extend over part of the model only. In order to 
model continuous layers, these dropped layers may be modelled with zero thickness. A 
layered model showing dropped layer model with four layers, the second of which is dropped 
over part of the model, as shown in figure 1.1. 



Layer 2 
4 is dropped 



Figure 1.1. Layered model showing dropped layer 

For each layer, the following properties are specified in the element real constant table of 
ANSYS, (Main Menu> Pre-processor> Real Constants) (accessed with REAL attributes), 
o Material properties (via a material reference number MAT) 
o Layer orientation angle commands 
o Layer thickness 

I.4.I.I. Material Properties 

As with any other element, the material properties command (Main Menu> Pre¬ 
processor Material Props> Material Models > Structural > Linear > Elastio Isotropic or 
Orthotropic) is used to define the linear material properties, and the TB command is used to 
define the nonlinear material data tables The only difference is that the material attribute 
number for each layer of an element is specified in the element's real constant table. The 
linear material properties for each layer may be either isotropic or orthotropic. Typical fibre- 
reinforced composites contain orthotropic materials and these properties are most often 
supplied in the major Poisson's ratio form. Material property directions are parallel to the 
layer coordinate system, which is defined by the element coordinate system and the layer 
orientation angle and is described below. 



1.4.1.2. Layer Orientation Angle 

This defines the orientation of the layer coordinate system with respect to the element 
coordinate system. It is the angle (in degrees) between X-axes of the two systems. By default, 
the layer coordinate system is parallel to the element coordinate system. All elements have a 
default coordinate system which can be changed using the ANSYS element attribute (Main 
Menu> Pre-processor > Meshing > Mesh Attributes > Default Attributes). 

1.4.1.3. Layer Thickness 

If the layer thickness is constant, then specify the thickness at node. Otherwise, the 
thicknesses at the four corner nodes must be given as input. Dropped layers may be 
represented with zero thickness. 

1.4.2. Defining the Constitutive Matrices 

This is an alternative to specify the individual layer properties and is available as an 
option [KEYOPT (2)] for SOLID46 and SHELL99. The matrices, which represent the force- 
moment and strain-curvature relationships for the element, must be calculated outside the 
ANSYS program. 

The main advantages of the matrix approach are: 

1. It allows you to incorporate an aggregate composite material behaviour. 

2. A thermal load vector may be supplied. 

3. The matrices may represent an unlimited number of layers. 

The terms of the matrices are defined as real constants. Mass effects are incorporated 
by specifying an average density for the element. If the matrix approach is used, detailed 
results in each layer cannot be obtained since individual layer information is not input. 



1.5. Sandwich and Multiple-Layered Structures 

Sandwich structures have two thin faceplates and a thick which is relatively weak core and 
the same is illustrated in Figure 1.2 . 

Faceplate 

Core (at least 1/2 of total thickness) 

Faceplate 



Figure 1.2. Sandwich Construction 

The sandwich structures can be modelled with SHELL63, SHELL91, or SHELL181; 
SHELL63 has one layer but permits sandwich modelling through the use of real constants. 

1.5.1. Node Offset 

For SHELL91 and SHELL99 the node offset option locates the element nodes at the 
bottom, middle or top surface of the shell. The figures 1.3 shows layered shell with nodes at 
mid plane and Figure 1.4 shows the layered shell with nodes at bottom surface, and these 
surfaces are aligned. 



Figure. 1.3. Fayered shell with nodes at mid plane 






Nodes located on bottom surface witfi KEYOPTfl 1) = 1 


Figure 1.4. Layered Shell with Nodes at Bottom Surface 

1.5.2. Specifying Failure Criteria 

Failure criteria are used to learn if a layer has failed due to the applied loads. The 
failure criteria can be chosen from three predefined criteria. 

The three predefined criteria are: 

1. Maximum Strain Failure Criterion, which allows nine failure strains. 

2. Maximum Stress Failure Criterion, which allows nine failure stresses. 

3. Tsai-Wu Failure Criterion, which allows nine failure stresses and three additional 
coupling coefficients. 

1.6. Introduction to Composites 

Composite materials are formed by combining two or more materials that have quite 
different properties. The different materials work together to give the composite unique 
properties, but within the composite, the different materials are apart and they do not dissolve 
or blend into each other. Technically, the composite has both good compressive strength and 
good tensile strength . 

The threads of glass in fibreglass are very strong under tension but they are also brittle 
and will snap if bent sharply. The matrix not only holds the fibres together, it also protects 
them from damage by sharing any stress among them. The matrix is soft enough to be shaped 
with tools, and can be softened by suitable solvents to allow repairs to be made. Any 
deformation of a sheet of fibreglass necessarily stretches some of the glass fibres, and they 
are able to resist this, so even a thin sheet is very strong. It is also quite light, which is an 
advantage in many applications. 
















Over recent decades many new composites have been developed, some with very 
valuable properties. By carefully choosing the reinforcement, the matrix, and the 
manufacturing process that brings them together, engineers can tailor the properties to meet 
specific requirements. They can also select properties such as resistance to heat, chemicals, 
and weathering by choosing an appropriate matrix material. 

1.6.1Selection of material for Matrix 

For the matrix, many modern composites use thermosetting or thermo softening 
plastics (also called resins). The use of plastics in the matrix explains the name 'reinforced 
plastics' commonly given to composites. The plastics are polymers that hold the 
reinforcement together and help to determine the physical properties of the end product. 
Thermosetting plastics are hquid when prepared but harden and become rigid (ie, they cure) 
when they are heated. The setting process is irreversible, so that these materials do not 
become soft under high temperatures. These plastics also resist wear and attack by chemicals 
making them very durable, even when exposed to extreme environments. Thermo softening 
plastics, as the name implies, are hard at low temperatures but soften when they are heated. 
Ceramics, carbon and metals are used as the matrix for some highly specialised purposes. For 
example, ceramics are used when the material is going to be exposed to high temperatures in 
heat exchangers and carbon is used for products that are exposed to friction and wear in 
bearings and gears. 

1.6.2 Materials for the reinforcement 

Although glass fibres are by far the most common reinforcement, many advanced 
composites now use fine fibres of pure carbon. Carbon fibres are much stronger than glass 
fibres, but are also more expensive to produce. Carbon fibre composites are light as well as 
strong. They are used in aircraft structures and in sporting goods (such as golf clubs), and 
increasingly are used instead of metals to repair or replace damaged bones. Even stronger 
(and more costly) than carbon fibres are threads of boron. Polymers are not only used for the 
matrix, they also make a good reinforcement material in composites. For example, Kevlar is a 
polymer fibre that is immensely strong and adds toughness to a composite. It is used as the 
reinforcement in composite products that require lightweight and reliable construction i.e., in 
structural body parts of an aircraft. Composite materials were not the original use for Kevlar 



- it was developed to replace steel in radial tyres and is now used in bullet proof vests and 
helmets. 

1.6.3. Manufacturing process 

Making an object from a composite material usually involves some form of mould. 
The reinforcing material is first placed in the mould and then semi-liquid matrix material is 
sprayed or pumped in to form the object. Pressure may be applied to force out any air 
bubbles, and the mould is then heated to make the matrix set solid. The moulding process is 
often done by hand, but automatic processing by machines is becoming more common. One 
of the new methods is called pultrusion (a term derived from the words 'puli' and 'extrusion'). 
This process is ideal for manufacturing products that are straight and have a constant cross 
section, such as bridge beams. 

In many thin structures with complex shapes, such as curved panels, the composite 
structure is built up by applying sheets of woven fibre reinforcement, saturated with the 
plastic matrix material, over an appropriately shaped base mould. When the panel has been 
built to an appropriate thickness, the matrix material is then cured. 

In many advanced composites , the structure may consist of a honeycomb of plastic 
sandwiched between two skins of carbon-fibre reinforced composite material. Such sandwich 
composites combine high strength, and particularly bending stiffness, with low weight and 
costly. 

1.6.4. Advantages of composites 

• The greatest advantage of composite materials is strength and stiffness combined with 
lightness. By choosing an appropriate combination of reinforcement and matrix 
material, manufacturers can produce properties that exactly fit the requirements for a 
particular structure for a particular purpose. 



• Composites are less likely to break up completely under stress than metals (such as 
aluminium). A small crack in a piece of metal can spread very rapidly with very 
serious consequences (especially in the case of aircraft). The fibres in a composite act 
to block the widening of any small crack and to share the stress around. 

• The right composites also stand up well to heat and corrosion. 

• Another advantage of composite materials is that they provide design flexibility. 
Example: a surfboard or a boat hull. 

1.7. Classification of materials 


Degree of 

anisotropy 

Principal 

axes 

Properties 

Example 

Isotropic 

Orthogonal 

Constant regardless of 

direction 

Metals 

Square symmetric 

Orthogonal 

Two different 

principal axes 

Unidirectional fibres or 

woven cloth 

Orthotropic 

Orthogonal 

Three different 

principal axes 

Unidirectional weave with 

light weft 

Anisotropic 

Any angle 

Constant relative to 

axes 

Filament wound tube : 

Many crystals 

Aeolotropic 

Any angle 

May change with 

position 

Timber 





















1.8. Properties of composites 

The different properties of the composites are listed below. 

1.8.1. Fibre volume fraction (Vf): 

It is calculated by using the formula 



Pjt 


where: 

n = the number of layers, 

Af = the areal weight of the fabric, 
pf= density of the fibre, and 
t = the thickness of the laminate 

1.8.2. Areal Weight of a Fabric (Af) 

It is calculated by formula 

Ay = 2 N f N T 7TfP f 

where, for a balanced fabric, the parameters are: 

Nf = number of filaments per tow 
Nt = number of tows in unit width of fabric 
rf = radius of the fibre cross-section 
pf= density of the fibre 

1.8.3. Rule of Mixtures 

The properties of composite can be calculated from rule of mixtures. It is given by 

E ' = w, v f E t +v A 

where 

E c = Young's modulus of the composite 
Ef = Young's modulus of the fibre 
E m = Young's modulus of the matrix 
Vf = fibre volume fraction 
V m = matrix volume fraction (1-Vf-V v ) 

V v = void volume fraction 



r|i = fibre length distribution factor 
r| 0 = fibre orientation distribution factor 


1.8.4. Stacking sequence 

There are a number of ways in which fibres can be arranged in order of increasing 
stiffness and strength. 

• 3-D random: e.g. injection moulding grades. 

• planar random: e.g. moulding compounds, chopped strand mat and 
continuous random swirl. 

• quasi-isotropic: e.g. continuous fibres oriented at 0/-45/90/+45 or 0/60/120. 

• bidirectional: e.g. woven fabrics or cross-plied unidirectional at 0/90. 

• unidirectional: e.g. pultrusions and aligned fibre laminates. 

The normal way to concisely record a laminate stacking sequence is as follows: 

[07+457-45790°] ns 

where the subscripts are: 

n is the number of times the sequence is repeated 
s means the laminate is symmetric 
T is the total number of plies, and 

overbar denotes that the laminate is symmetric about the mid-plane of the ply 
Thus for n = 2 , when * denotes the line of symmetry, the sequence will be: 
07+457-45790707+457-45790° * 907-457+45707907-457+4570° 




1.8.5. Young's modulus 

Hooke's Law originally stated that the stress (a) is proportional to the strain (s) with 
the constant of proportionality known as Young's modulus (E), so that a = Es. Traditional 
structural engineering materials are isotropic and hence the value of E is independent of 
direction. However, the Young's modulus of a composite material is anisotropic (varies with 
direction) and can be estimated using the rule-of-mixtures. The variation of young’s modulus 
for fibre orientation is as shown in figure 1.1 


s c = ViV^s B s + VA 


where: 


E c = Young's modulus of the composite 
Ef = Young's modulus of the fibre 
E m = Young's modulus of the matrix 
Vf = fibre volume fraction 
V m = matrix volume fraction (1-Vf-V v ) 

V v = void volume fraction 

rji = fibre length distribution factor 

i)„ = fibre orientation distribution factor 

1.8.6. The fibre length distribution factor (FLDF) 

It assumes that the matrix and fibre remain elastic and the interface bond is perfect, 
the shear stress at the fibre ends is maximum and falls to zero after half the "critical length", 
and the tensile stress at the fibre ends is zero and rises to a maximum after half the "critical 
length". 

FLDF can be calculated using equation : 



tanh ( jSL / 2) 





E f A f In [rJRj,) 


where 




G m = the shear modulus of the matrix, 

L = the fibre length, 

Af = the cross sectional area of the fibre, 
Rf = the radius of the fibre and 
R = the mean separation of the fibres. 


1.8.7. Critical length 

The critical length (when there is no debonding) is given by the expression 


} _ Rf °11 

I 

^12 


where R = the fibre radius, 

an = the tensile stress in the fibre and 

a '12 = the shear strength (of the interface or of the matrix as appropriate). 

1.8.8. Fibre orientation distribution factor 

The fibre orientation distribution factor can be calculated 

180 " 

= cos4 3 

i-0" 

For a unidirectional ply, the rule-of-mixtures does assume that the fibres are continuous with 
uniform cross-section and lie parallel to each other. There is perfect bonding between the 
fibre and the matrix. 
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Fig.1.1. Variation of young’s modulus with fibre orientation. 


1.8.9. Shear modulus 

The inter-relationship of the elastic constants for isotropic materials is given by 
G = E/2(l+v) 


K= E/3(l-2v) 

where E = the Young’s modulus, 

G = the shear modulus, 

K= the bulk modulus and 
v = Poisson’s ratio. 

For orthotropic materials within a single plane, the shear modulus is 

^ _ yjEA 

-\/p12 M21) 


1.8.10. Poisson's ratio 


Poisson's ratio is denoted by the Greek letter v. It has a value determined by: 
v = -(strain normal to the applied strain)/(strain parallel to the applied strain). 














Tensile deformation is taken as positive and compressive deformation is taken as negative. 
The minus sign in the definition of Poisson's ratio normally results in positive values of 
Poisson's ratio and a thermodynamic constraint which restricts the values to -1 < v < 1/2. 

1.8.11. Bulk modulus 

The bulk modulus for a square symmetric material can be found by the following 
equation 

g _ 1-^2-^3 

3(l - l^ 1 l^ 3 ) 

This equation reduces to the isotropic form when Ei = E 2 = E 3 and V12 = V21 = V23. 

1.9. Research problem and objectives: 

Besides the applications in the aeronautical and automobile industry, the application 
of composite materials has been enlarged, including new areas such as nautical industry, 
sporting goods, civil and aerospace construction. In order to find the right combination of 
material properties and in service performance, their dynamic behaviour of the composites 
need to be studied. To avoid structural damages caused by undesirable vibrations, it is 
important to determine the natural frequencies of the structure to avoid resonance, the mode 
shapes to reinforce the most flexible points and the right positions to reduce weight. With 
respect to these dynamic aspects, the composite materials represent an excellent possibility to 
design beams with requirements of dynamic behaviour. 





1.9.1. Objectives: 


• The main objective of this thesis is to contribute for a better understanding of 
the dynamic behaviour of doubly symmetric composite I- beams. 

• To know the influence of fibre orientation, elastic modulii and height to 
thickness ratio on natural frequencies of composite I- beams for uniform and 
tapered beams for various boundary conditions. 

• To investigate the dynamic behaviour of the beams, experimental studies and 
numeric analysis using Finite element method have been carried out and the 
results are presented in detail. The experimental results of a cantilever glass- 
epoxy rectangular composite beam are used to validate the results obtained 
from finite element analysis . 

1.10. Scope of the Current Research Work 

Since many years, the combination of different materials has been used to achieve 
better performance requirements. Although the benefits brought by the composite materials 
known , the present focuses the behaviour of composite I-beams is dealt in detail, and the 
influence of fibre orientation, elastic modulii and height to thickness ratio on natural 
frequencies of composite beams. The modal analysis and harmonic analysis is done by 
changing the fibre angle in top and bottom flanges, the effect of fibre angle rotation in web of 
simply supported composite I-beams. The effect of Elastic modulus ratio and height to 
thickness ratio of composite beams on the natural frequencies for a simply supported and 
clamped composite beams were investigated with the stacking sequence of the top and 
bottom flanges at [0/90] keeping the web unidirectional. 

Further, the study was extended to know the effects of fibre orientation, modulus ratio 
for a cantilever composite I - beam with a small end fixed on natural frequencies and mode 
shapes. The natural frequencies and mode shapes of the same cantilever I beam were 
obtained when its bigger end is fixed, with the same parameters. The dynamic analysis is 
extended to know the behaviour and mode shapes of a composite I-beam with circular holes 
in web with the same dimensions, parameters and boundary conditions. 



The flow charts 1.8 shows the numerical analysis and experimental analysis plan of 
verifying the results and validation for rectangular glass-epoxy composite laminates. The 
flow chart 1.9 shows the study of orthotropic I-beams by different methods namely Galerkin 
Finite Element Analysis, Dynamic stiffness method and Finite element analysis using 
ANSYS with various boundary conditions. The flow chart 1.10 shows the detailed plan of 
total research work in thesis. 



1.11. Flow chart showing the numerical and experimental plan 














































1.13. Flow chart showing the detailed plan of thesis 













































1.14. Flow chart to show the Procedure to do Modal analysis using 
ANSYS 

Choose the Element type 

4 

Choose the Appropriate Real Constants and Material Constants 
Create key points at the start and end of the beam 
Create a line between the above points 

Mesh the beam with the size control 

Apply Boundary Conditions 

Change the analysis type to MODAL 

Solution > Analysis Type >New Analysis >Modal 

G 

Set the Analysis Options 
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Solution > Analysis options >Extract 10 mode>OK 

Enter <50> for the ending frequency 

Solve the load Set 

G 
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1 

List RESULTS summary 

General Post Processor >List Results>Results Summary 


Read Results for a Substep 



General Post Processor > Read results > First Set 
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PLOT deformed geometry 

ft 

Read Results for Next Substep 

General Post Processor > Read results > Next Set 



1.15. Flow chart to show the Procedure to do Harmonic analysis of a 
Cantilever beam using ANSYS 


Choose the Element type 

Choose the Appropriate Real Constants and Material Constants 
Create key points at the start and end of the beam 
Create a line between the above points 

Mesh the beam with the size control 

Apply Boundary Conditions 
Apply 1.0 N Load at the free end 
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New analysis > Harmonic 
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Set the Solution Method 
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V 

Set the Tolerance to 1.0 e -8 < OK> 


V 
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i 

Solution > Load step Options - Time Frequency > Frequency and Sub Steps 


V 



Set the Harmonic Frequency Range 0-50 Hz 


V 

Set the Number of Sub steps to 100 

V 

Set to be Stepped 

V 

SOLVE the model 



1.16. Organisation of the thesis 

The thesis is organized in the following manner. 

Chapter-1 covers the importance, background and the present status of research in the fields 
relevant to the present work. It starts with the importance of thin walled isotropic beams, 
assumptions made to develop a model for composite beams in particular for I-section. The 
flow charts clearly explain the steps adopted for experimental analysis , orthotropic I-beams 
by different methods, presentation of work in thesis, modal and harmonic analysis of 
cantilever beams. 

Chapter-2 presents the literature survey on finite element modeling and simulations and 
vibrations of composite beams of open and closed section including isotropic and composite 
beams. This chapter reviews the extent of literature that is relevant to the problem, which the 
author intends to investigate. About 131 papers were referred and are given in reference 
section at the end of the thesis. In particular, the literature pertaining to composite beams of 
open and closed sections was considered. 

Chapter-3 covers the scope of the research work which includes modelling of isotropic 
beams resting on elastic foundation. Basic governing equations for various boundary 
conditions are derived in this chapter. The mode shapes for Indian standard I- beams 
(uniform and tapered) are modelled using ANSYS and are shown. 

Chapter 4 illustrates the experimental set up to conduct impact test for cantilever rectangular 
glass-epoxy composite laminated beam .The experiment is conducted with glass-epoxy 
specimens with different stacking and results are obtained and tabulated. These results are 
validated by finding the percentage variation with numerical results obtained by ANSYS 
.Therefore the numerical values are in good agreement with experimental results. Hence, the 
research work was extended for composite I-beams (both uniform and tapered beams) 
assuming the above validation. 

Chapter 5 emphasises Galerkin’s finite element formulation of orthotropic I-beams and 
governing equations for tapered cantilever isotropic I-beam. Analytical results are presented 
for various different conditions. 



Chapter 6 presents a dynamic stiffness matrix method has been developed for computing the 
natural frequencies and buckling loads of thin-walled beams resting on continuous elastic 
foundation and can be applied for treating beams non-uniform cross section and also non- 
classical boundary conditions. 

Chapter 7 presents mathematical modeling of composite I-Beam and the Stiffness Equation 
of a K th orthotropic lamina in the laminate coordinate system are developed and explicit 
expressions for I-section are given. Finite element formulation and the closed form 
expressions for vibration frequencies of a simply supported anisotropic beam, expressions for 
axial frequency, flexural frequencies in x and y directions, and the torsional vibration 
frequency are presented. 

Chapter 8 covers the effect of various parameters namely fibre orientation, modular ratio and 
height to thickness ratio on natural frequencies of composite I- beams and mode shapes are 
shown for cantilever boundary condition. The harmonic analysis is done for doubly 
symmetric composite I-beam. 

Chapter 9 deals with modelling of tapered composite I-beam with its small end fixed and 
modal and harmonic analysis are carried out using ANSYS R10.0, the analysis software. The 
results obtained are presented in graphical form. 



Chapter 10 deals with modelling of tapered composite I-beam with its small end free and big 
end fixed .Modal and harmonic analysis are carried out using ANSYS R10.0, the analysis 
software. In this chapter, a case study is done for a cantilever composite I-beam with three 
holes of diameter 32 mm in web. Modal and harmonic analysis is done and compared. 

At the end of the thesis, discussion on results, conclusions, future scope of work, 
author’s index, and references are presented. 
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Chapter 2 

LITERATURE REVIEW 

In the literature survey, it is found that thin-walled beams of open cross-section, 
such as I-section beams, are susceptible to instability in a variety of modes, but a few 
publication have dealt with dynamic behaviour of such members. Closed-form solution for 
flexural and torsional natural frequencies of isotropic thin-walled beams is found in the 
literature [4-6]. However, the research was limited to doubly symmetric beams with simply 
supported boundary conditions. Free vibration of thin-walled beams with variable I-section 
by finite element method was studied by Wekezer [4,5] . Ohga et al. [7] presented an 
analytical procedure to estimate the natural frequencies and mode shapes of thin-walled 
members using the transfer matrix method and various types of cross-sections of thin-walled 
members including C-section, box-section and I-section were analyzed. 

The literature mentioned in the list dealt with free vibration of isotropic metallic thin- 
walled beams. For composite thin-walled beams, the flexural and torsional vibrations are 
fully coupled in general, even for a doubly symmetric cross-section due to their material 
anisotropy. Only a few works have addressed the dynamic behaviour of composite thin- 
walled members. Bauld and Tzeng [8] extended Vlasov's thin-walled bar theory [9] to 
symmetric fibre-reinforced laminates to derive buckling equations of composite thin-walled 
beams. In addition , the differential equations for the general case of tapered wide -flanged 
beam columns have been derived using the Vlasov’smethod[9] for uniform beams. An 
experimental investigation of the elastic stability of tapered beam columns has also been 
investigated. Song and Librescu [10] proposed an analytical model to predict free vibration 
behavior of laminated composite thin-walled beams of open sections. They investigated 
natural frequency and mode shape with respect to the fibre orientation for composite box- 
section beams. Recently Lee and Kim [11] presented a general analytical model to predict 
buckling loads and mode shapes of thin-walled composites. 

Xi-Xian Wu and C.T.Sun[13] gave the simplified theory for composite beams in 1992. 
Torsional vibrations of beams and cracked beams were dealt in detail by J.M.Gere[14]. 
Pratyoosh Gupta, Wang S.T, Blandford, G.E[15] gave the detailed study of Torsional Buckling 
of Non-Prismatic I-Beamand derived the equations for various boundary conditions. 



Shakourzadeh ,,et al [16] deduced torsion bending element for thin-walled beams with open 
and closed cross sections. He also studied the torsional vibration of beams with doubly 
symmetric cross section for which the shear centre and centroid coincide and analysed the 
effect of warping on the frequencies of torsional vibrations and shapes of the normal modes 
of vibration for bars of single span with various conditions. Derivation of the fundamental 
theory of strength and stability of thin-walled members was performed Xie (20) presented 
nonlinear normal modes of a clamped -clamped beam for large amplitude displacements. 
Naidu ,chen,Banerjee [21-24] have studied the free vibration of thin walled beams with axial 
loading and extended this to beams subjected to torsion and bending. Banerjee [25] has 
developed the dynamic stiffness matrix method for an axially loaded Timoshenko beam 
element. He has considered the case of coupled bending and torsion for the beam in his basic 
governing differential equations. The detailed procedure for analysis of orthotropic and 
anisotropic beams using various elements listed in ANSYS [26] library. 

The results obtained for the fundamental normal mode are compared with the 
corresponding reported experimental and theoretical studies. Godier, Timoshenko,Vlasov and 
others. Free Torsional Vibrations of Tapered Cantilever I-beams and its formulation is given 
by Kameswara rao and Mirza [38].Timoshenko[52] initiated the concept of non uniform 
torsion and considered warping of cross sections of a symmetrical I-beam subjected to torsional 
moment. Wagner [53] generalised the Timoshenko torsion theory. The finite element method 
has been shown by Barsoum [54], and Barsoum and Gallangher [62] ,to be completely 
general, in that they had provided results that are sufficiently accurate for engineering 
purposes. 

Gere and Carter[57] obtained the critical buckling loads for tapered columns. A finite 
element formulation using Galerkin’s method for the buckling problem of tapered members 
was presented by Morrel and Lee. Gere and Lin [56] generalised the theory of vibrations of 
thin-walled beams of an arbitrary open section. Goodier [58] published a series of studies in 
which he simplified and proved some of the assumptions proposed by earlier investigations. 
Nowisinki[60] and Panovko[61] have studied several literature papers related to the theories of 
thin-walled beams in which they have focussed on optimum design of tapered thin-walled 
beams. 


H.N.Lee [63] presented an analysis of non-uniform torsion of tapered I-beam in 1956, 
the taper being only of a restrictive type. Using the finite element technique, Krajcinovic[ 64] , 



developed a formulation for thin-walled members based on the hyperbolic functions to 
express the twist. These functions which are the solutions to the exact differential equation 
for twist, lead to complicated stiffness expressions in torsional and warping constants. It does 
not include the effects of instabilities due to torque. The formulation is only applicable to 
solid beams where the shear centre coincides with the centre of gravity. The investigations 
done prior to the year 1945 related to lateral stability and flexural -torsional stability of uniform 
thin-walled beams were compiled and consolidated by Timoshenko [68]. Valsov’s extensive 
investigations of thin-walled elastic members were published in book form in 1940. A new 
edition containing comprehensive study of equilibrium, stability and vibrations of thin-walled 
members of arbitrary cross sections was published in 1958. 

Static torsional response and lateral-torsional stability of tapered I-beams has been 
investigated by many researchers [65-72]. Hamaychi [66], Lee and Szabo [67], Wilde [68] 
presented basic derivations for a comprehensive theory of non-uniform torsion of tapered I- 
beams. Among these, Massey and McGuire [70] studied the lateral stability of stepped 
cantilever beams of rectangular and I-cross-section using a Runge Kutta integration 
procedure. The problem of lateral-torsional buckling of tapered I-beams has been studied by 
Kitipomchai and Trahair [71], Brown [72], Culver and Preg [46] and Shiomi and Kurata [47]. 
In these studies, solutions were obtained for simply supported and cantilever beams using 
finite-difference or finite-integral methods.A review of the literature clearly shows that very 
few studies [73-74, 38] have been conducted on the free vibration characteristics of non- 
uniform thin-walled beams. Several investigations [84-86] have been reported on torsional 
vibrations and stability of long uniform thin-walled open section beams, such as I-beams. 

Wiedemann (76) was presented generalised method and analytical solution is also 
presented for the natural frequencies ,mode shapes and orthogonality conditions of an 
arbitrary system of Euler-Bernowli’s beams inter connected by arbitrary joints and subjected 
to boundary conditions. Prokic (77) derived the governing equations for coupled bending and 
torsional vibrations for simply supported beams. Claeyssen (79) determined the modes of 
uniform beam described by Euler-Bernoulli model with axial force or Timoshenko Model in 
a unified manner with the use of a fundamental free response. Several authors have studied 
vibration of beams masses. Goel (80) investigated the free vibration of a beam -mass system 
hinged at either end by rotational springs and arbitrarily located heavy mass. Chang (82) 



analysed the vibrations of a mass loaded system with a heavy tip mass by using Laplace 
Transform. 

Lin (83) derived the exact governing equations for linear vibrations of a rotating 
Timoshenko beam by the d’Alemberts principle and Virtual Work Principle, on the 
assumptions that the beam is linear elastic and the steady state axial strain is small. 
Prokic(84)considered the problem of thin walled beams of closed sections by means of exact 
solution and the differential equations are derived from the principle of virtual work due to 
variation of displacements, for the case of simply supported beam. Gurgose (84) studied the 
natural frequencies of restrained beam and rods with point masses. Dynamic stiffness matrix 
of a moving beam is developed and investigated the free flexural vibration characteristics by 
Banarjee(85), 

Law, S.S., (86) investigated the dynamic behaviour of a bolted joint which has 
flexibility in both the tangential & rotational directions the joint is prestressed with a axial 
tension in the volt shank. The formulation of a hybrid beam column element including the 
end springs is presented, and the dynamic behaviour of a cantilevers beam with nonlinear 
semi rigid joint is studied, The natural frequencies & node shapes remain relatively 
unchanged only for a limited range of the joint stiffness. Mei, C. (88) presented the wave 
vibration analysis of axially loaded bending torsion coupled composite beam structures, 
including the effects of axial force, shear deformation and rotary inertia & this study was 
limited for a axially loaded Timoshenko composite beam. The study also includes the 
material coupling between the bending tensional modes of deformation that is usually present 
in laminated composite beam due to ply orientation. 

The classical single-mode elliptic function frequency solutions for simply supported 
beams and square plates are assessed and the characteristics of non linear response are 
studied. Magadi Mohareb,M(89) , developed a finite element formulation that captures both 
the St.Venant and warping torsional effects of open sections. The formulation is applicable to 
prismatic thin-walled open sections commonly used in steel construction and is able to 
represent general torsional and warping end restraints. Gere,J.M.(90) studied the fee torsional 
vibrations of bars of thin walled open cross section for which the shear centre and centroid 
coincide, for I-beams and Z-sections. The differential equations for torsional vibrations are 
derived including the effect of warping. The effects of warping on frequency of vibration are 




determined for bars of single span with various end conditions. For a simply supported bar, a 
formula for the principal torsional frequencies and an expression for mode shape are derived. 

Prasad, K.S.R.K.,(91) studied the analysis of Eigen value problems such as free 
vibration problems, and exact solutions of the governing differential equations are possible 
only in special cases and presented a comparative study of some commonly used approximate 
methods. A comparative study of Raleigh Ritz method and iterative type approximate 
methods for the evaluation of Eigen values has been made with the aid of examples. Kwok- 
Tung Chen (92) derived the equation of free transverse vibration of beams with two sections 
of partially distributed mass and its exact solution has been obtained. Experimental data for a 
cantilever beam are given to verify the computational results. Yuchengshi(93), presented a 
multimode time-domain modal formulation based on the finite element method for large 
amplitude free vibrations of thin composite plates. Banerjee (94) presented and exact stiffness 
matrix for a composite Beamed including the effects of shear deformation and rotatory 
inertia, for a Timoshenko beam. An explicit analytical expression for each of the elements of 
the dynamic stiffness matrix is derived by use of REDUCE computing package. Hamed, E 
(95) investigated analytically dynamic behaviour of reinforced concrete (RC) beams 
strengthened with externally bonded composite materials. The equations of notion along with 
the boundary & continuity conditions are derived using Hamilton’s variational Principle and 
the Kinematic relations of small deformations the mathematical formulation also includes the 
constitutive laws that are based on beam & laminations theories, and the two - dimensional 
elasticity representation of the adhesive layer including the closed form solution of its 
stresses displacement. 

Banerjee[95] has developed the exact dynamic stiffness matrix method of approach for 
composite Timoshenko beams and has considered all boundary conditions. Arpaci,A(96) 
presented an exact analytical method for predicting natural frequencies of beams with thin- 
walled open cross sections having no axis of symmetry. The governing equations give a 
characteristic equation of the 12 th order with real coefficients. The roots are found 
numerically and the exact boundary conditions are considered especially for free ends to 
obtain natural frequencies. Li [97] has studied the free torsional vibrations of open cross 
section and also derived expressions including the effect of warping. Abbas[98] has extended 
the study of vibrations of Timoshenko beams with elastically restrained ends and Lee S Y[99] 
extended his study for non uniform beams with elastically restrained ends. Bishop and Price 



(100), and Banerjee and Williams (101) have taken into the effects of rotary inertia and 
transverse shear. Warping stiffness has been neglected in both works. A recent study carried 
on by Bercin and Tanaka (102) has included the effect of warping in doubly coupled 
vibration analysis of Timoshenko beams with monosymmetric open cross sections. In case of 
triple coupling, where the flexural vibrations in two mutually perpendicular directions and the 
torsional vibrations are all coupled is dealt with by Arpaci and Bozdag(103). This is 
experienced when the cross sections of no symmetry are of concern. 

Kameswara Rao [104-112] has studied the torsional vibrations and stability of lengthy 
thin walled beams on elastic foundation. He has developed the governing equations for the 
short wide flanged thin walled beams of open sections. He also studied the effects of 
longitudinal inertia and shear deformation. He has derived the equations to find vibration 
frequencies and mode shapes for generally restrained Bernouli beams. Then, he deduced the 
equations for frequency analysis of Clamped-Clamped uniform beams with intermediate elastic 
support. The same concept was extended to two span uniform Euler-Bernouli beams. He has 
also considered non linear torsional vibration of thin walled beams of open sections. 

To avoid structural damage caused by undesirable vibrations and resonance, it is 
important to determine the natural frequencies of the structure. The mode shapes must also be 
determined in order to reinforce the most flexible points or to determine the right positions 
where is necessary to reduce weight or to increase damping [113]. Sahu and Asha (2005) 
[114], used eight-noded isoparametric quadratic shell element to develop the finite element 
procedure, to investigate the vibration and stability behaviour of pretwisted panels, the effect 
of various geometrical parameters like angle of twist, aspect ratio, lamination parameters, 
shallowness ratio are studied. Tita (2003) [115] worked on the theoretical and experimental 
dynamic analysis of reinforced epoxy resin. He used laminates in his study. The laminates 
were fabricated by hand-lay-up process and cut to beam shape specimens. The specimens 
were used as free-end beams in the vibration measurements. He calculated the mechanical 
properties of the composite analytically and used them in his simulations of the dynamic 
properties. He presented his ANSYS® simulation results in contour format showing the mode 
types and shapes. 

Colakoglu (2006) [116] performed vibration experiments on 10-layer beam specimens 
of glass polyethylene composite at a range of temperatures. The vibration was induced by the 
impact of spherical steel ball hammer. ANSYS® numerical simulations were also used to 



obtain the frequency response. Teng and Hu (2001) [117] analyzed the design parameters for 
constrained layer damping structures by employing the Ross-Kerwin-Ungar (RKU) model. 
They also discussed the effects of temperature, frequency and the dimensions of damped 
structures on vibration damping characteristics. 
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Chapter 3 

VIBRATIONS OF A THIN-WALLED I-BEAM 
RESTING ON ELASTIC FOUNDATION 

3.1. Introduction 

Static and dynamic analysis of beams resting on elastic foundation occupies a 
prominent place in contemporary structural mechanics. The vibrations and buckling of 
continuously supported finite and infinite beams resting on elastic foundation has an 
application in the design of highway pavements, aircraft runways and in the use of metal rails 
for rail road tracks. The vibrations of beams resting on elastic foundation have been 
investigated by a number of investigators. In these studies, it has been mostly assumed that 
the foundation reacts in tension as well as compression. Various types of foundation models 
such as Winkler, Pasternak, Vlasov, Filonenko-Borodich, etc., have been used in the analysis 
of structures on elastic foundation. Among these , the Winkler model, in which the medium is 
taken into account as a system composed of infinitely close linear springs, is the simplest and 
is adopted for analysis. It assumes that the foundation applies only a reaction normal to the 
beam that is proportional to the beam deflection. 

3.2. Derivation of basic differential equation: 

As the cross sectional dimensions are assumed to be small compared to the length in 
thin-walled beams, the second order effects such as longitudinal inertia and shear 
deformation can be treated as negligible. 

In this work, based on Timoshenko’s torsion theory, the governing differential 
equation of free motion of a doubly symmetric thin-walled I-beam on elastic foundation 
subjected to a time- invariant axial compressive load is derived based on Hamilton’s 
principle. 

According to Hamilton’s principle, 

( Tk - U + W )dt = 0 (3.1) 

Where Tk = kinetic energy of the strained bar, 

U = total strain energy, 

W = the potential energy of the external force, and 
to,t, are fixed constants. 

According to Saint-Venant, the cross- sections are assumed to rotate about the 



(3.2a) 


Centroid -shear centre ’O’, giving rise to a torsional couple, 

T s = G C s (&|>/Sz) 

Where G = shear modulus, 

Cs = torsion constant for the cross section, 
cp(z,t) = the angle of twist 
The torsion constant for an I- section is 

C s = (2bt f 3 + ht w )/3 (3.2b) 

Where b= width of the flanges, 
t, =thickness of flanges 

h = height between the centre lines of the flanges. 
t w = thickness of the web. 

The strain energy U 1 at any instant t in the beam of length L due to saint - Venant torsion is 
Ui = (Vi) J GC S (d<\>/dz ) 2 dz (3.2c) 

Accompanying the rotation in a warping of the section which is assumed constant in each 
piece of cross- section having a moment M. and the X- displacement of the top flange 
center line, U, is given by, U = (h/2) (j) (3.2d) 

And hence the moment M in the top flange is given by 

M = E If (d 2 u/dz 2 ) = E If (h/2)(o 2 <j)/Sz 2 ) (3.2e) 

Where E = Young’s modulus, 

I f = moment of inertia of each flange area about y-axis. 

The shear force Q due to the bending of the flanges is given by 

Q = (dM/dz) = E If (h/2)(aV0z 3 ) (3.2f) 

The equal and opposite shear forces Q, a distance h apart in the top and bottom flanges, give 
rise to torque due to warping, T w , is given by 

T w = - Qh = - E If (h 2 /2)( d 3 <\>/dz 3 )= - E C w ( d 3 $/dz 3 ) 
where C w = If (h 2 /2) = warping constant for an I- section 


(3.2g) 



The total torque, T t , on the cross section is given by 

T t = T s + T w = GCs( dtydz) - E C w ( d^ld z 3 ) (3.2h) 

If U 2 is the strain energy of the two flanges due to warping, then 

U 2 = Vi j 2EI f ( S 2 (j)/Sz 2 ) 2 dz = Vi J E C w ( d 2 tyd z 2 ) 2 dz (3.2i) 

The strain energy U 3 due to the Winkler type elastic foundation, is given by 

U 3 = 1 / 2 J K t ((j)) 2 dz (3.2j) 

Hence the total strain energy U, at any instant t becomes 

u = Ui + U 2 + U 3 = Vi I GCs ( dfy/dz) 2 + E C w ( 5 2 ^/0z 2 ) 2 + K t (c^) 2 |dz (3.2) 

The kinetic energy of rotation of the cross section at the corresponding instant t is given as: 
T = I /2 p I p (d§/dt) 2 dz (3.3) 


Where I p is the polar moment of inertia of the cross section and 
p the mass density of the material of the beam. 

The potential energy due to the external time - invariant axial compressive load, P, acting at 
the centroid of the cross section at the corresponding instant t is given by 

W = Vi J (Pip /A) (dtydz f dz (3.4) 

Where A is the area of the cross section 

Substituting for Tk, U and W from the above equations in the equation (1), taking 
thevariations of the integrand and integrating the first term by parts with respect to ‘t’ and the 
next four terms with respect to ‘z’, we get, 



(3.5) 



Assuming that the values of 0 given at the two fixed instants, the second integral 
vanishes. If the boundary conditions are such that the third and fourth integrals also vanish, 
then the associate differential equation of motion is given by: 



- Kt(</>) - pip 
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(3.6) 


3.3. Natural boundary conditions 


In deriving the basic differential equation of motion (6) from equation.(5)it was assumed that 
the expressions 
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(3.8) 


Equations (3.7) and (3.88) give the natural boundary conditions for the finite bar, and are 
satisfied if the end conditions are taken as 


(l)^ = Oand ( ^ 


dz‘ 


O 


(3.9) 


These conditions imply restraint against rotation but not against warping; that is, the 
end of the bar does not rotate but is free to warp. This is the case of a “simple support”. 



(2) (j) =0 and (dtydz) = 0 


(3.10) 


These conditions imply restraint not only against rotation but also against any warping of the 
end cross section this means that the end of the bar is built in rigidly so that no deformation 
of the end cross section can take place. These conditions define a “fixed support”. 



these conditions imply no restraint of any kind at the end of the bar. This requires that the 
bending moment in the flange ends and torque acting on the end cross-section must be zero. 
These correspond to a “free end”. 



The later conditions imply no warping and zero shear forces in the end flanges. These 
conditions are useful for finding symmetric modes of vibration in simply supported, fixed- 
fixed and free- free beams. 

3.4. Time- dependent boundary conditions 

The homogeneous boundary conditions discussed above, give the free vibrations of 
bars. For forced vibrations produced by the motion of boundaries, approximate time 
dependent end conditions are given by prescribing at each end one member of each of the 
Products: 



Or equivalent of M5 (dtyd z) and T t Sc)) respectively. 

Of the many conditions thus obtained, the following are of more theoretical interest: 

1. Twisting moment T t prescribed, flange bending moment M = 0 or (S<j)/Sz) = 0. 

2. (j) or (3(j)/5t) prescribed, flange bending moment M = 0 or (S<j)/5z) = 0. 



3. Flange bending moment M prescribed, twisting moment T t =0, cj)=0, 

(dfy/dz) or (d 2 $/dzdt) prescribed, twisting moment T t =0 or (j)= 0. 

In the case of semi-finite beams, conditions need be prescribed at one end since all physical 
quantities at any instant are zero at the free end. 

3.5. Analysis of various terms 


1. If K t = P = 0 and C w = 0 , equation (6) becomes 
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This equation represents Saint Venant’s torsion theory for slender beams and does not 
include warping of the cross section shear deformation and or longitudinal inertia effects. 


2. IF K t = P = 0 The EQN.(6) reduces to 
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(3.14) 
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This equation represents Timoshenko’s torsion theory which includes the effect of warping of 
the cross section. 


3. If K t = 0, Equation.(6) reduces to 
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(3.15) 


This Equation represents the effect of an axial time- invariant compressive load added to 
Timoshenko’s torsion theory. 



4. If P = 0, Equation.(6) reduces to 

Gcl d V 2 1 - EcJ d ^/ 4 ) - - pIpi^A . 2 1 = 0 (3.16) 

\ /dz ) \ /dz J V /dt 

This equation represents the effect of Winkler type constant modulus elastic foundation 
added to Timoshenko’s torsion theory. 

3.6. Non-Dimensionalization and General Solution of Equation of 
Motion 

For mathematical simplification, it is convenient to reduce equation (6) to a non- 
dimensional form, simultaneously introducing some dimensionless parameters having 
physical interpretations. 

Introducing, Z = z/L, the non-dimensional beam length, and 

Ti = (EC w )/(pIpL 4 ) t, the dimensionless time variable, 

Equation (3.16) in non-dimensionless form can written as: 

a V 4 V (^ 2 - V 2 / 5 V 2 ] + V (</>) + f= 0 (3.17) 

\ /dz ) l /dz J V /dt J 

Where k 2 = (GCL 2 )/(EC W ) , warping rigidity parameter (3.18) 

V 2 = (PI p L 2 )/(AECw) , axial load parameter (3.19) 

and y 2 = (K t L 4 )/(4EC W ), foundation parameter, (3.20) 

the general solution of equation.(3.14) can be obtained by using the standard method of 

separation of variables. Thus, by taking (j) in the form 

(j) = X (Z) T (ti) (3.21) 

and then substituting into equation.(3.14),separating the variables, and setting the resulting 

expressions equal to Zn 2 , we obtain 

T = A„ cos A n ti + B„ Sin A n ti (3.22) 

The expression for a normal mode of vibration is then 

(j) = X (A„ cos A„ ti + Bn Sin A„ ti) (3.23) 

in which X is the normal function giving the shape of the mode of vibration. 


An is the dimensionless torsional frequency parameter given by 
A n 2 - (pIpL 4 p n 2 )/(ECw) 


(3.24) 



A,, is the dimensionless torsional frequency parameter. Any actual motion of the vibrating 
beam can be obtained by a summation of normal modes, so that in the general case 
(f>\ 1 = X \X n (A n cos A rl ti + B n Sin A„ ti) 

(3.25) 

in which the coefficients A n and B n are found from the initial conditions of the vibration. 
The equation for determining the normal function X, found by substituting equation (3.24) 
into differential equation(3.14), then 

( a, y 8z M k2 ~ v2 P 2 y 8z M 4 r i -^ 2 )x=o ( 3 . 26 ) 

The general solution of this equation may be found by taking the normal function X in the 
form: 


X = D e nz 

The equation (3.27) yields the auxiliary algebraic equation. 


(3.27) 


if - (k 2 - V 2 )(ij 2 )+ (4y 2 - An 2 ) = 0 (3.28) 

The four roots of the equation are 

7i = +«,; v 2 = -« 2 ; % = +Waba = (3-29) 

In whichai and Pi are the positive, real quantities given by 

“i = ( 1 / • jy {(£ 2 - V 2 ) + (£ 2 - V 2 ) 2 + 4 ( A „ 2 - 4/ 2 )] 1,2 } 1/2 (330) 
A = (1/72)“{ -(k 2 ~V 2 ) + (k 2 -V 2 ) 2 +4U„ 2 - V] 1 ' 2 )} 1 ' 2 (331) 

The general solution of equation (26) then becomes either 

X = D y e a ' z + D 2 e~ a ' z + D 2 e p ' z + D 4 e~ filZ 

or 

X — D x cosh a x z + D 2 sinh a x z + D 3 cos J3 x z + D^Sinp x z ( 3 . 32 ) 


There are four arbitrary constants in this expression, which must be determined so as to 
satisfy the particular boundary conditions of the problem. For any beam there will be two 
boundary conditions at each end and these four conditions determine the frequency equation 
and the ratios of three of the constants to the fourth constant. Solving the frequency equation 



then determines the principal frequency of vibration. With the frequencies and normal 
functions determined, the solution is essentially complete. 

3.7. Frequency Equations and Modal Functions 

Because of the complexity of the frequency equations, the study is limited the case of 
simply supported beam. 

3.7.1. Boundary conditions 

In terms of non-dimensional parameters, the boundary conditions can be written as: 

1. Simple support: 

X = 0 , d 2 X/dZ 2 = 0, (3.33) 

2. Fixed support: 

X = 0, dX/dZ = 0, (3.34) 

3. Free end: 

d 2 X/dZ 2 = 0 , (k 2 - V 2 )dx/dz - d 3 x/dz 3 = 0 (3.35) 

3.8. Simply Supported Beam: 

The present study is limited to simply supported beams only. The boundary 
conditions are applied and frequency equations and mode shapes for simply supported beam 
are established in detail. From the boundary conditions 

X = d 2 X/dZ 2 = 0 at Z=0, and X = d 2 X/dZ 2 = 0 at Z = 1. 

For the conditions at z=0, equation (32) gives: 

D 3 + Di = 0, 
and Di(ar +P 1 2 ) = 0 . 

Since the (ar +P 1 2 ) is not equal to zero, it follows that D 3 = Di= 0. 

From the second pair of conditions, equation (32) gives: 

D 2 Sinh ai + D 4 Sin Pi =0, and D 2 0.2 Sinh ai - D 4 P 2 Sin Pi =0, 

The characteristic equation is (ar + Pr)Sinh ai sin Pi =0 

The frequency equation is Sin Pi = 0, Therefore Pi = nn , where n=l, 2, 3... 

The frequency parameter is 



A = 


n 


1—2 / 2_2 

7i {n n 


+ k' 


■a 2 )+V 


1/2 


3.9. Finite Element Simulation 


Finite Element modelling and analysis is used to analyze natural frequencies and 
mode shapes of Indian Standard Junior Beam (ISJB - 150), Indian Standard Light Beam 
(ISLB - 150) and tapered beams with following dimensions by using Finite Element 
Analysis. 


Problem description of an ISJB -150 Beam 

Weight per meter, W = 7.1 kg-f 
Sectional Area, a = 9.01 cm 2 
Depth of Section, h = 150 mm 
Width of Flange bf= 50 mm 
Thickness of Flange, tf = 4.6 mm 
Thickness of Web, t w = 3.0 mm 
Young’s Modulus, E = 2E5 
Poison’s Ratio = 0.3 

Problem description of an ISLB -150 Beam: 

Weight per meter, W = 14.2 kg-f 
Sectional Area, a = 18.08 cm2 
Depth of Section, h = 150 mm 
Width of Flange, bf = 80 mm 
Thickness of Flange, tf = 6.8 mm 
Thickness of Web, t w = 4.8 mm 
Young’s Modulus, E = 2E5 


Poison’s Ratio = 0.3 



3.9.1 Meshing model of an ISJB-150. 




Fig.3.1 Meshing model of an ISJB-150 

3.9.I.I. Different mode shapes of an ISJB-150. 
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Fig.3.2 Mode shapes of an ISJB-150 
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Table.3.1. furnishes the frequencies of first 10 mode shapes 
of an ISJB-150. 
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Fig. 3.3. Frequency Variation of an ISJB -150 Beam Over Different Set 








3.9.2. Meshing model of an ISLB-150. 
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Table 3.3. The frequencies of first 10 mode shapes of ISJB-150 
tapered beam. 


SET 

FREQUENCY, 

Hz 

1 

0.14934E-01 

2 

0.28969 E-01 

3 

0.42584 E-01 

4 

0.68405 E-01 

5 

0.10413 

6 

0.11395 

7 

0.14493 

8 

0.16825 

9 

0.16956 

10 

0.23494 
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FREQUENCY 


2.50E-01 




Fig. 3.9 shows frequency variation of an ISJB - 150 tapered beam. 










3.9.4. Meshing model of an ISLB-150 tapered beam. 
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Fig.3.10 Meshing model of an ISLB-150 tapered beam. 
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3.9.4.I. Different mode shapes of an ISLB-150 tapered beam: 
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Fig.3.11 Different mode shapes of an ISLB-150 tapered beam 
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Table 3.4. The frequencies of first 10 mode shapes of ISLB-150 tapered 
beam. 


SET 

FREQUENCY, Hz 

1 

0.23549E-01 

2 

0.40699 E-01 

3 

0.60454 E-01 

4 

0.95891 E-01 

5 

0.14304 

6 

0.14833 

7 

0.18048 

8 

0.20608 

9 

0.20863 

10 

0.21864 
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Figure. 3.12shows frequency variation of an ISLB - 150 tapered beam 
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SETS 


Fig. 3.13 shows the comparison of variation of frequency over different sets of all four beams. 
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3.10. RESULTS AND DISCUSSIONS 

It has found that the frequencies of standard beams are higher than that of tapered beams, the frequencies of individual beams 
gradually increases over different sets and also numerical results obtained from the ANSYS software are almost nearer to the 
theoretical values obtained by using wittrick.W.H. & Williams, E.W., Algorithm. 
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3.9.2.I. Different mode shapes of an ISLB-150: 
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Fifth mode shape 
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Tenth mode shape 

Fig.3.5 Mode shapes of an ISLB-150 




























Table 3.2. The frequencies of first 10 mode shapes of ISLB-150 


SET 

TIME/FREQ 

1 

0.14501E-01 

2 

0.29164 E-01 

3 

0.93341 E-01 

4 

0.13815 

5 

0.16644 

6 

0.20584 

7 

0.22566 

8 

0.23144 

9 

0.24260 

10 

0.24351 
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Figure 3.6. Frequency variation of an ISLB -150 Beam over different sets 












3.9.3. Meshing model of an ISJB-150 tapered beam: 
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Fig.3.7 Meshing model of an ISJB-150 tapered beam 
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3.9.3.I. Different mode shapes of an ISJB-150 tapered beam: 
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Chapter 4 

EXPERIMENTAL ANALYSIS OF GLASS-EPOXY COMPOSITE 
RECTANGULAR LAMINATED BEAMS 


4.1. Introduction 

It is known from dynamic systems that a structure temporarily excited by an external force will vibrate at a frequency defined as the 
natural frequency. The natural frequency is characteristic of the entire structure and its boundary conditions. An external excitation applied 
to the structure at the same frequency as the structure's natural frequency will result in resonance. At resonance, the structure will vibrate at 
higher than normal amplitude levels. Depending upon the overall system design, the amplitude of vibration at resonance can cause unwanted 
operating conditions as well as a catastrophic failure. This type of system failure can be avoided by determining the structure's natural 
frequency and maintaining the overall system at an operating condition which will not induce resonance. The natural frequencies of a 
structure can be determined by measuring frequency response. One means of measuring the frequency response is to impact test the 
structure. Impact testing is a measurement of the ratio of the beam's transient response to the excitation impulse versus frequency. 

4.2. Manufacturing of Glass-Epoxy Composites: 

Manufacturing of composite materials involves distinct operations that may vary depending upon available technology, existing 
facilities and skilled persons. The manufacturing process may also vary due to wide variety of composite materials and their 
applications. Factors considered for selection of most efficient manufacturing process are 
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• User needs 

• Performance requirements 

• Labour 

• Surface complexity 

• Volume of the production 

• Size of the product 

• Appearance 

• Production rate 

• Tooling/assembly 

• Equipment 

The main goals of the composite manufacturing process are 

• To achieve a consistent product by controlling fibre thickness, fibre volume and fibre direction. 

• To minimise voids. 

• Reduce internal residual stresses. 

• Economy of process 

The hand lay-up method is used for manufacturing the glass - epoxy composite specimens for conducting the experiments. Even 
though the method has been replaced with automated techniques, the lay-up of preimpregnated material by hand is the oldest and most 
common fabrication method for advanced composite structures. The following steps illustrate the adopted manufacturing process for 
glass epoxy specimens. 
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1. The surface of the tool is cleaned and release agent is applied. The release agent can be in liquid or solid form. 

2. An optical sacrificial layer is laid upon the tool surface. This layer is usually a fibre glass fibre made with the same resin 
system as the composite laminate. The sacrificial layer protects the laminate from surface abrasion and surface irregularities 
during manufacturing. 

3. A peel ply is placed on top of the sacrificial layer. The peel ply is removed after processing. 

4. The preimprignated plies are cut according to design specifications. They can be cut by hand using shears or a steel blade 
knife. However , automated cutting machines have largely replaced hand cutting. 

5. The first prepeg ply is oriented and placed upon the tool or mold. Subsequent plies are placed one upon another; a roller is used 
to compact the plies and remove entrapped air that could later lead to voids. 

6. A flexible resin dam is anchored to the sacrificial layer approximately 3mm from the edge of the laminate. The dam prevents 
resin flow out of the laminate, in the plane of the laminate. 

7. Another peel ply is placed on top of the laminate to protect the laminate surface. 

8. A sheet of porous release film is laid over the dam and the laminate. The porous release film will serve as a barrier to prevent 
bonding of the composite laminate to the secondary materials to follow. 

9. Next, bleeder plies are laid up over the release film. The bleeder plies extend to the edge of the laminate. 

10. Another porous release ply is next laid up over the bleeder plies. 

11. Caul plates are sometimes placed on the top of lay-up. The caul plate is steel or aluminium plate that protects the surface from 
sharp temperature increases (it acts as a heat sink) and it gives a smooth non-waxy surface texture. 

Obviously, hand lay-up represents an extreme that can not be scaled to high volume or large components. There is a significant 
amount of skilled labour necessary for the hand lay-up of composite parts. Each step has a specific purpose and function. This type of 
fabrication is the most time-consuming, but it is also the most flexible manufacturing process. 
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4.3 Dimensions of the rectangular glass -epoxy composite beam 

Two different symmetrical laminates with twenty layers each is manufactures by laying one layer over another by pressure. 
Type 1: [45/-/45/45/-45/45/-45/45/-45/0/90] s 
Type 2: [0/90/0/90/0/90/0/90/0/90] s 

The details of the Rectangular beam considered for both experimental and simulation by ANSYS R10.0 are : 

Beam length = 40 cm , 

Width of the beam = 25 mm, 

Thickness of the beam = 1.6mm, and 
Total mass of the beam = 28xl0" 3 kg. 

The rectangular cantilever beam is modelled and analysed using ANSYS R 10.0 with the above data using SHELL99, for both the 
types of stacking of layers. 

Material used: glass fibre and Epoxy resin (HY -956) 

Table 4.1.Material properties 


Material 

Properties 

Symbol 

Value 

Glass fibre 

Elastic modulus 

E f 

76xl0 9 


Density 

pr 

2.56xl0 3 
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Poisson’s ratio 

V12 

0.22 

Epoxy resin 

Elastic modulus 

En 

4.0xl0 9 


Density 

pm 

1 . 3 x 103 


Poisson’s ratio 

Vm 

0.40 

Laminate 
{orthotropic) 

Elastic modulus 

Ell 

44.8xl0 9 

Fibre direction 

Normal to fibre 

E 22 =E33 

11.27xl0 9 


Density 

Pc 

1780 


Shear modulus 

G12= G13 

4.86xl0 9 



G23 

4.45xl0 9 


Poisson’s ratio 

V12 

0.28 


Fibre volume fraction 

Vf 
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4.4. Modal Analysis by Finite Element method 

The glass -epoxy composite rectangular beam with above dimensions is modelled using SHELL99. This beam id discretized 
into 50 elements and analysis is done by ANSYS software by fixing one end. The analysis is carried out for two different symmetrical 
laminates with twenty layers each by laying one layer over another with following orientations. 

Type 1: [45/-/45/45/-45/45/-45/45/-45/0/90] s 
Type 2: [0/90/0/90/0/90/0/90/0/90] s 

4.5. Mode shapes: 

Few of the mode shapes obtained from Finite Element Analysis are shown below for both the cases. 

4.5.1. Mode shapes for first case: 
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Case 1 —Fourth mode 
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4.6. Experimental analysis of Glass-Epoxy composite Rectangular 
laminated beams: 

The experiment is conducted by fixing the specimen in a rigid support with other end free to vibrate, as a cantilever beam. The 
impact hammer was used to give impact to the specimen, and the Spectral analyzer was set from 0-400 Hz. This output was captured 
by using accelerometer connected to the sensor which is fixed to the beam and amplified .The Frequency response function (FRF) is 
analysed using the spectrum analyzer. After the measurement of the FRF (Amplitude and phase), the natural frequencies were 
evaluated. The block diagram of the experimental set up is as shown in figure 4.1. 
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1. Cantilever beam 
4. Power Amplifier 

Figure.4.1. 


2. Accelerometer 3. Impact hammer with piezoelectric cells 
5. Spectrum Analyzer. 

Block diagram of experimental set up for modal analysis 
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4.7. Background 

Frequency response measurements are a type of modal testing which can be used to determine the resonant frequencies of a 
structure. Frequency response testing is accomplished by exciting a system with a known input and simultaneously measuring the 
corresponding output response. Figure 4.2 is a block diagram which represents a system with one input (single-point excitation) and 
the output response. 


X(tJ) 



► Y(u) 


Fig. 4.2.Block diagram with input and output response. 


Where X(ro) is the excitation supplied to the structure, Y(cd) is the system's response, and H(co) is the physical system being 
analyzed. The ratio of response to excitation versus frequency is defined as the frequency response.Figure4.3 shows the instruments 
that can be used to conduct frequency response measurements by impact test of a beam and utilizing digital data acquisition and set up 
is shown. The fig.4.4. shows the photograph of two specimens of Glass-Epoxy rectangular composite specimens, which are used for 
conducting the experiments to find the natural frequencies. These specimens are manufactured by following the steps mentioned in 
this chapter. 
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Fig.4.3.Picture of the experimental set up of an impact experiment. 
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Fig.4.4. photograph of specimens of Glass-Epoxy rectangular composite specimens. 







4.8. Comparisons of experimental results with numerical results 

The rectangular laminated beam is analysed by ANSYS for a cantilever boundary condition and numerical results are listed in 
table 4.2 for both the cases. The rectangular laminated beam is fixed at one end in the impact experimental setup and results are 
tabulated in table 4.3 for both cases of beams. 

Table: 4.2Numerical results obtained from FEM 


Case 1 

[45/-45/45/-45/45/-45/45/-45/0/90] 

Case 2 

0/90/0/90/0/90/0/90/0790 

Mode 

Frequency in Hz 

Mode 

Frequency in Hz 

1 st Flexural Mode 

4.2 

1 st Flexural Mode 

5.3 

2 nd Flexural Mode 

26.2 

2 nd Flexural Mode 

34.1 

3 rd Flexural Mode 

74 

3 rd Flexural Mode 

96.5 
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Table: 4.3.Experimental results for two cases 


Case 1 

[45/-45/45/-45/45/-45/45/-45/0/90] 

Case 2 

0/90/0/90/0/90/0/90/0790 

Mode 

Frequency in Hz 

Mode 

Frequency in Hz 

1 st Flexural Mode 

4 

1 st Flexural Mode 

5.9 

2 nd Flexural Mode 

27 

2 nd Flexural Mode 

36 

3 rd Flexural Mode 

74 

3 rd Flexural Mode 

97 

4 th Flexural Mode 

149 

4 th Flexural Mode 

109 

1 st Torsional 

Mode 

181 

1 st Torsional 

Mode 

197 

5 th Flexural Mode 

237 

5 th Flexural Mode 

326 


The natural frequencies obtained from numerical analysis and experimental are compared separately and the percentage 
variation is calculated and listed in table 4.4 and table 4.5 for case 1 and case 2 of rectangular laminated composite beams. This is to 
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validate the results obtained from numerical analysis. It is found that the natural frequencies obtained from ANSYS are in good 
agreement with the experimental results. The same is assumed to be valid for the cases of uniform composite I-beam and tapered 
composite I-beam and the study of these beams is carried out in subsequent chapters. 
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Table: 4.4.Percentage Variation in natural frequencies for case 1 


Case 1 

[45/-45/45/-45/45/-45/45/-45/0/90] 

Difference 

% 

Variation 

Mode 

Frequency in 

Hz from 

FEM 

Frequency in 

Hz from 

experiment 

1 st Flexural 

Mode 

4.2 

4 

0.2 

4.7619 

2 nd Flexural 

Mode 

26.2 

27 

0.8 

3.053 

3 rd Flexural 

Mode 

74 

74 

0.0 

0.0 

4 th Flexural 

Mode 

143 

149 

6.0 

4.196 
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Table: 4.5.Percentage Variation in natural frequencies for case 2 






4.9. Results and Discussions 


The natural frequency values for the composite glass-epoxy rectangular laminated beams are shown in the above tables 4.2 and 
4.3 .These results are compared with the numerical results obtained from ANSYS and the percentage variation is found to be the in 
good agreement with the experimental values of frequencies. 
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Chapter 5 

Galerkin Finite Element Analysis of Free Torsional Vibrations 
of Tapered Cantilever I - Beams 
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Chapter 5 


Galerkin Finite Element Analysis of Free Torsional Vibrations of Tapered 

Cantilever I - Beams 


5.1. Introduction 

It is commonly believed that material saving can be accomplished by using non-uniform beams. Non-uniform thin-walled 
beams are quite commonly used in aircraft, bridges and several other industrial structures. In view of the current emphasis on 
structural optimization, there is an urgent need to obtain a proper understanding of the vibration characteristics of non-prismatic 
structural members. Several investigations [84-86] have been reported on torsional vibrations and stability of long uniform thin- 
walled open section beams, such as I-beams. 

Static torsional response and lateral-torsional stability of tapered I-beams has been investigated by many researchers [65-72]. 
Among these, Hamaychi [66], Lee [67], Wilde [68]], presented basic derivations for a comprehensive theory of non-uniform torsion of 
tapered I-beams and they have also studied the lateral stability of stepped cantilever beams of rectangular and I-cross-section using a 
Runge Kutta integration procedure. The problem of lateral-torsional buckling of tapered I-beams has been studied by Kitipornchai and 
Trahair [71], Brown [72], Culver and Preg [46] and Shiomi and Kurata [47]. In these studies, solutions were obtained for simply 
supported and cantilever beams using finite-difference or finite-integral methods. 
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A review of the literature clearly shows that very few studies [73-74, 75] have been conducted on the free vibration 
characteristics of non-uniform thin-walled beams. It can be also seen that especially little progress has been made in the area of 
torsional vibrations of I-beams with taper along their depth. 

In this chapter, the governing differential equation for torsional vibrations of tapered doubly symmetric I-beams is derived and 
solved for the case of a cantilever I-beam fixed at its bigger end and free at the larger end by utilizing the Galerkin finite element 
method. The governing differential equation for torsional vibrations of tapered doubly symmetric I-beams was derived and solved and 
analysed for cantilever I-beam. The individual as well as combined effects of linear taper in the width of flanges and the depth of the 
web on torsional natural frequencies were investigated. 

This chapter deals with the case of cantilever I-beams with uniform linear taper in web height ((3), flange width (a) and flange 
thickness (y) with the bigger end fixed and smaller end free. 


5.2. Governing equation of torsional motion 

Consider free torsional vibrations of a doubly symmetric thin-walled I-beam of variable cross-section and of length L. the 
variations in the width of the flanges and the depth of the web are respectively assumed to be of the form 

b(z) = b 0 (l -a/L), (5.1) 

d(z) = do (1 - pz/L), (5.2) 
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where z is the distance along the length of the beam, and a and p are the taper ratios in the width of the flanges and depth of the web, 
respectively. The values with subscript zero are those at the bigger end of the beam . 

When the beam is executing torsional vibration, the bending moment, M, induced in the flanges is given by equation 5.3. 

M (z) =EI f (z)0 2 u/5z 2 . (5.3) 

In equation (5.3), the warping displacement, u, in the flanges is given by 

U(z) = [h(z)/2](K (5.4) 

Where h(z) is the height between centrelines of the flanges. Hence, 

H(z) = d(z) + t f . (5.5) 

Substituting equation (5.4) in equation (5.3) gives 

M(z) = Elf (z) e 2 /5z 2 [(h (z)/2) (5 2 u / dz 2 )] (5.6) 
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The torque, T, induced in the beam is given by (5.6) 


T(z) = Gc s (z) 0(() / dz - h(z) 3M (z)/ dz + M(z) dh(z) / dz. (5.7) 

For free torsional vibrations, the static torque is replaced by the inertia torque which has an intensity of pl p (z) o 2 (f)/ot 2 , 
the polar moment of inertia and p is the mass density of the material of the beam. Thus, 

5T (z)/ dz = pl p (z) d 2 $/dt 2 . (5.8) 

Substituting equation (5.7) in equation (5.8) gives 

d/dz [Gc s (z) d§ldz\- <3 2 /<3z 2 [ EIf(z)h(z) <3 2 /0z 2 ((h(z)/2)((|))] 

+2(0/0z)[EIf(z) (dh(z)/dz)/dz (d 2 /dz 2 (h(z)/2) ()>)] = pI P (z) (0 2 (j)/5t 2 ) (5.9) 

For harmonic vibration the angle of twist (j) can be written in the form 

(j) (z,t) = 0 (z) sin cot. (5.10) 


Substituting equation (5.10) in equation (5.4) gives 



did z [Ecw(z) d 2 /dz 2 (h(z)/2) 0]- d/dz[Gc s (z) d0/dz ] 

-2(d/dz)[EIf(z) (dh(z)/dz)/dz (d 2 /dz 2 (h(z)/2) 0)] - pl p (z) co 2 0 =0 (5.11) 

When expanding the differential equation (5.11), terms containing d 2 h(z)/dz 2 are neglected as the depth taper considered in this study 
is linear. 

Equation (5.11) is a fourth order differential equation with variable coefficients and a Galerkin-finite element method has been utilized 
in this study in obtaining the solutions for cantilever beams. 

5.3. Finite Element Formulation 

In the Galerkin-finite element formulation, the domain of the beam is subdivided into a number of elements. A third degree 
polynomial for angle of twist in z is assumed over the element. The variation of Torsional amplitudes 0 e over each element, in terms 
of nodal degrees of freedom 0 and 0', is given by 

{ 0e} — [x] { §e} (5.12) 

Where {x} is the row matrix of shape functions 
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(5.13) 


{5} e = [01, 0*1, 02, 0'2] T 


Here primes denote differentiation with respect to z and subscripts 1 and denote the two ends of the element. 
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Substituting equation (5.12) in equation (5.11) gives the residual R e for the element as 


R e = Elf h [h/2] " [ X] " {5} e - [Gc s X] {5} e - 2 EI f h' [h/2] " [ X] " {5} e -pI P co 2 x{5} e 

(5.14) 

Where E = Young’s modulus, 

If = Moment of Inertia of flange 
h = height of the beam, 

[h/2] = second derivative of h w.r.t z 
[ X] = second derivative of X w.r.t z 
h' = first derivative of h w.r.t z 

In the Galerkin - finite element method, the weighted residual is minimized by setting 

a L 

_1 R e 0edz = O (5.15) d{5}e 0 

This procedure is repeated for all elements and after the usual assembly procedure; the final matrix equation governing free vibration 
is obtained as 
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[A]{ 5} -X 2 n [B]{ 5} = 0. 


(5.16) 


where 


A 2 n = Afoco 2 L 4 /EI f0 (5.17) 

the quantities Afo and lfo are the flange area and moment of inertia respectively at the fixed end. Also, [A] is the stiffness matrix, [B] 
the mass matrix and {5} the eigenvector. Equation (5.16) can be solved by using any standard eigenvalue algorithm to obtain the 
eigenvalues and eigenvectors {5}. 


5.4. Numerical Results 

Numerical results of the frequency parameter A, have been obtained, by dividing the beam into a sufficient number of elements, for the 
case of a tapered cantilever I-beam. It is found that 12 elements are sufficient to yield convergence. The boundary conditions 
considered at the fixed end are 0 = 0' =0. 

The values of A, for various combinations of taper ratios a and p for the first three modes of vibration are presented in the figures 
Graph 5.1., Graph 5.2.and Graph 5.3 (first, second and third modes respectively). The combined influence of positive taper ratio P and 
negative taper ratio a on the first three frequencies of vibration are shown in the figures. 
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ALPHA 


Fig. 5.1.Combined influence of positive taper ratio p and negative ratio 

(FIRST MODE) 
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Fig. 5.2 Combined influence of positive taper ratio p and negative ratio a 
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(SECOND MODE) 


Fig. 5.3Combined influence of positive taper ratio p and negative ratio 

(THIRD MODE) 



5.5. Results and discussion 


The torsional frequencies of first three modes of vibration are increasing significantly by increasing values of taper ratio [1 
and marginally for increasing values of taper ratios a and A,. 

The torsional frequencies of first three modes of vibration are significantly increasing for increasing values of taper ratios [1 
and X and are decreasing drastically for increasing values of taper ratio a in second mode and third modes of vibration only. 

The torsional frequency increases significantly for increasing values of combined influence of negative taper ratios a and p in 
first mode. In second and third modes of vibrations the torsional natural frequency gradually decreases to certain value and then starts 
increasing. 
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Chapter 6 


Torsional Vibrations and Buckling of Thin-Walled Beams on Elastic Foundation- 

Dynamic Stiffness Method 

6.1. Introduction 

The problem of free torsional vibration and buckling of doubly symmetric thin- walled beams of open section, subjected to an 
axial compressive static load and resting on continuous elastic foundation, is discussed in this chapter. An analytical method based on 
the dynamic stiffness matrix approach is developed including the effect of warping. The resulting transcendental equation is solved for 
thin-walled beams clamped at one end and simply supported at the other. The dynamic stiffness matrix can be used to compute the 
natural frequencies and mode shapes of either a single beam with various end conditions or an assembly of beams. When several 
elements are to be used the over all dynamic stiffness matrix of the complete structure must be assembled. The associated natural 
frequencies and mode shapes are extracted using Wittrick-Williams algorithm (36). The algorithm guarantees that no natural 
frequency and its associated mode shape are missed. This is, of course, not possible in the conventional finite element method 
Numerical results for natural frequencies and buckling load for various values of warping and elastic foundation parameter are 
obtained and presented. 

The vibrations and buckling of continuously - supported finite and infinite beams on elastic foundation have applications in 
the design of aircraft structures, base frames for rotating machinery, rail road tracks, and etc. several studies have been conducted on 
this topic, and valuable practical methods for the analysis of beams on elastic foundation have been suggested. A discussion of 
various foundation models was presented by Kerr. 
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While there are a number of publications on flexural vibrations of rectangular beams and plates on elastic foundation, the 
literature on torsional vibrations and buckling of beams on elastic foundation is rather limited. Free torsional vibrations and stability of 
doubly-symmetric long thin-walled beams of open section were investigated by Gere , Krishna Murthy and Joga Rao and Christiano 
and Salmela , Kameswara Rao et al., used a finite element method to study the problem of torsional vibration of long thin-walled 
beams of open section resting on elastic foundations. In another publication Kameswara Rao and Appala satyam developed 
approximate expressions for torsional frequency and buckling loads for thin walled beams resting on Winkler-type elastic foundation 
and subjected to a time invariant axial compressive force. 

It is known that higher mode frequencies predicted by approximate methods arte generally in considerable error. In order to 
improve the situation, a large number of elements or terms in the series have to be included in the computations to get values with 
acceptable accuracy. In view of the same, more and more effort is being put into developing frequency dependent ‘exact’ finite 
elements or dynamic stiffness and mass matrices. In the present chapter, an improved analytical method based on the dynamic 
stiffness matrix approach is developed including the effects of Winkler - type elastic foundation and warping torsion. The resulting 
transcendental frequency equation is solved for a beam clamped at one end and simply supported at the other. Numerical results for 
torsional natural frequencies and buckling loads for some typical values of warping and foundation parameters are presented. The 
approach presented in this chapter is quite general and can be utilised in analyzing continuous thin - walled beams also. 

6.2. Formulation and analysis: 

Consider a long doubly-symmetric thin walled beam of open section of length L and resting on a Winkler -type elastic 
foundation of torsional stiffness Ks. The beam is subjected to a constant static axial compressive force P and is undergoing free 
torsional vibrations. The corresponding differential equation of motion can be written as 
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6.1 


ECw 0 4 (j)/0z 4 — (GCs -pIp/A)0 2 (j)/Sz 2 + K s ())+ pIp)o 2 (^/5t 2 

In which E, the modulus of elasticity; Cw ,the warping constant; G, the shear modulus; Cs , the torsion constant; Ip, the polar moment 
of Inertia; A, the area of cross section; p, the mass density of the material of the beam; (j),the angle of twist; Z, the distance along the 
length of the beam and t, the time. 


For the torsional vibrations, the angle of twist c()(z,t) can be expressed in the form 

c()(z,t) = x(z)e ipt 6.2 

in which x(z) is the modal shape function corresponding to each beam torsion natural frequency p. 

The expression for x(z) which satisfies equation(6.1) can be written as: 


x(z) = A Cos az + B Sinaz +C Coshpz +D sinpz 6.3 

in which 

aL.pL = (1/V2) { + (k 2 -A 2 ) + [(k 2 -A 2 ) 2 +4(Z 2 -4y 2 )] 1/2 } 1/2 6.4 

K 2 = L 2 GCs/ECw.A 2 =p Ip L 2 /AECw 6.5 

And A, 2 = pip L 4 p 2 n/ ECw y 2 =Ks L 4 /4ECw 6.6 
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From equation (4), the following relation between aL and pL is obtained. 

(PL) 2 = (cxL) 2 +K 2 -A 2 6.7 

Knowing a and p. The frequency parameter A, can be evaluated using the following equation: 

= (oL)(pL) + 4y 2 6.8 

The four arbitrary constants A, B, C, and D in equation (6.3) can be determined from the boundary equation of the beam. For 
any single span beam, there will be two boundary conditions at each end and these four conditions then determine the corresponding 
frequency and mode shape expressions. 

6.3. DYNAMIC STIFFNESS MATRIX: 

In order to proceed further, we must first introduce the following nomenclature: the variation of angle of twist (j) with respect to z is 
denoted by 0(z); the flange bending moment and the total twisting moment are given by M(z) and T(z). Considering clockwise 
rotations and moments to be positive, we have, 

0(z) = d(()/dz, 6.9 

hM(z) = -ECw( d2(|)/dz2) 
and 

T(z) = -ECw ( d>dz 3 ) + (GCs -pip/A) dcji/dz 6.10 
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Where ECw = I f W, 


If = the flange moment of inertia and 

h= the distance between the centre lines of flanges of a thin-walled I-beam. 

Consider a uniform thin-walled I-beam element of length L as shown in fig.l. 

By combining the equation (6.3) and (6.9), the end displacements (f)(0) and 0(0) and end forces, hM(0) and T(0) of the beam at z = 0, can 

be expressed as: 


^ <Kor 


r 

1 

0 

1 

0 


0(0) 

= 

0 

a 

0 

P 

B 

hM(0) 

>- < 

ECwa 2 
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-ECwp 2 

0 

C 

T(0) 

V. J 


0 

V 

ECwap 2 

0 

- ECwa 2 p 

D 
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Equation (6.11) can be written in an abbreviated form as follows: 

5(0) = V(0)U 6.11 


in a similar manner , the end displacements , (f(L)and 0(L) and end forces hM(L) and T(L), of the beam where z = 1 can be expressed 
as follows: 
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8(L) = V(L)U 
where 

{5(L)}T = {4>(L), 0(L), hM(L), T(L)} 


{U}T = {A,B,C,D} 
and 


r 


c 

-as 


s 

ac 


C 

ps 


s 

pc 


A 


[V(L)]=< 


ECwa 2 c ECwa 2 s -ECwp 2 C -ECwp 2 S 
- ECwap 2 s ECwap 2 c - ECwa 2 pS - ECwa 2 pC 




J 


in which c = CosaL; s = SinaL; C = CoshpL; S= SinhpL. 


By eliminating the integration constant vector U from equation (11) and (12), and designating the left end of the element as i and the 
right end as j. the equation relating the end forces and displacements can be written as: 


^ Ti 
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" jn 

J12 

J13 

J14 
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HMi 

j21 

j22 

j23 

j24 

0i 

Tj 

= j31 

j32 

j33 

j34 

<Pj 

HMj 

j41 

j42 

j43 

j44 

0j 


Symbolically it is written 

{F} = [J] {U} (6.12) 

where 

{F}T = {Ti, hMi, Tj, hMj} 

{U}T = {<|>i, 0j, c>j, 0j} 

In the above equations the matrix [J] is the ‘exact’ element dynamic stiffness matrix, which is also a square symmetric matrix. 

The elements of [J] are given by: 

jn = H[(a 2 +p 2 }(aCs+pSc] \ 
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ji 2 = -H[(a 2 -p 2 }(l-Cc)+2apSs] 


jn = - H[(a 2 +p 2 }(as+pS)] 
ji4 = -H[(a 2 +p 2 }(C-c)] 

j 22 = -(H/ap)[(a 2 +p 2 }(aSc-pCs)] -(6.13) 

j 2 4=(H/aP)[(a 2 +p 2 }(aS-ps)] 

j 23 = -j 14 
j 33 = j 11 
j 34 = “j 12 

j44 = j 22 

and H= EC w /[2ap[ (l-Cc)+( p 2 -a 2 ) Ss] 

using the element dynamic stiffness matrix defined by equation (6.12) and (6.13). One can easily set up the general equilibrium 
equations for multi-span thin-walled beams, adopting the usual finite element assembly methods. Introducing the boundary conditions, 
the final set of equations can be solved for eigen values by setting up the determinant of their matrix to zero. For convenience in 
programming, the signs of end forces and end displacements used in equation are all positive and are defined as shown in fig.6.1. 
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6.4 Method of Solution: 

Denoting the assembled and modified dynamic stiffness matrix as [DS], we state that det IDS I =0 

(6.14) 

Equation (6.14) yields the frequency equation of continuous thin-walled beams in torsion resting on continuous elastic 
foundation and subjected to a constant axial compressive force. It can be noted that equation (6.14) is highly transcendental in terms of 
eigen values X. The roots of the equation (6.14) can, therefore, be obtained by applying the Regula-Falsi method and the Wittrick - 
Williams algorithm on a high speed digital computer. Exact values of frequency parameter X for simply supported and built in thin- 
walled beams are obtained in this chapter using an error factor s= I O' 6 

6.5 RESULTS AND DISCUSSIONS: 

The approach developed in this paper can be applied to the calculation of natural torsional frequencies and mode shapes of 
multi -span doubly symmetric thin-walled beams of open section such as beams of I-section. Beams with non uniform cross sections 
also can be handled very easily as the present approach is almost similar to the finite element method of analysis but with exact 
displacement shape functions. All classical and non- classical (elastic restraints) boundary conditions can be incorporated in the 
present model without any difficulty. 

To demonstrate the effectiveness of the present approach, a single span thin walled I-beam clamped at one end (z=0) and 
simply supported at the other end(z=l)is chosen. The boundary conditions for this problem can be written as: 
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<K0) = 0; 0(0) = 0 


(6.15) 


(j)(l) = 0; M(l) = 0 (6.16) 

Considering a one element solution and applying the boundary conditions defined by equation (6.15) and (6.16) gives, 
j22=0 (6.17) 

This gives, 

j22 = -(H/a(3)[(a2 +(32}(aSc-(3Cs)] =0 (6.18) 

as H and (a2 +P2) are ,in general, non-zero, the frequency equation for the clamped, simply supported beam can , therefore, be 
written as 

atanhpL = ptanaL (6.19) 

equation (6.19)is solved for values of warping parameter k=l and k=10 and for various values of foundation parameter y in the range 
0-12. 

Figures 6.2 and 6.3 shows the variation of fundamental frequency and buckling load parameters with foundation parameter for 
values of k equal to 1 and 10 respectively. To give an idea about the accuracy of the results obtained, comparison is made with results 
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obtained using approximate methods such as Galerkin technique. It can be stated that even for the beams with non-uniform sections, 
multiple spans and complicated boundary conditions accurate estimates of natural frequencies can be obtained using the approach 
presented in this chapter. 

A close look at the results presented in figures clearly reveal that the effect of an increase in axial compressive load parameter 
Ais to drastically decrease the fundamental frequency A,(N=1). Further more, the limiting load where /.becomes zero is the buckling 
load of the beam for a specified value of warping parameter K and foundation parameter yone can easily read the values of buckling 
load parameter Acr from these figures for A,=0, as can be expected, the effect of elastic foundation is found to increase the frequency of 
vibration especially for the first few modes. However, this influence is seen be quite neglgible on the modes higher than the third. 
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Fig.6.1. First mode 
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Fig. 6.2 SECOND MODE 
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Buckling Parameter 
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Fig.6.4 FOUTH MODE 

Fig. 6. 5. Mode shapes for Wk= 0.01 
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Fig6.6. Mode shapes for Wk = 0.1 
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Fig.6. 7. Mode shapes for Wk = 1.0 
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6.6 CONCLUDING REMARKS: 


A dynamic stiffness matrix approach has been developed for computing the natural torsion frequencies and buckling loads of 
long, thin-walled beams of open section resting on continuous Winkler-Type elastic foundation and subjected to an axial time - 
invariant compressive load. The approach presented in this paper is quite general and can be applied for treating beams with non- 
uniform cross sections and also non-classical boundary conditions. Using Wittrick-Williams algorithm, a computer program has been 
developed to accurately determine the torsional buckling loads, frequencies and corresponding modal shapes .Results for a beam 
clamped at one end and simply supported at the other have been presented, showing the influence on elastic foundation, and 
compressive load and warping. While an increase in the values of elastic foundation parameter resulted in an increase in frequency, 
the effect of an increase in axial load parameter is found to be drastically decreasing the frequency to zero at the limit when the load 
equals the buckling load for the beam. 
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Chapter 7 

MATHEMATICAL MODELLING OF COMPOSITE I-BEAM 

7.1. Introduction to Composites 

The need for high performance to weight ratio structures coming from the most advanced engineering fields is the main 
driver of the increasing usage of composite materials for critical applications. Composite materials find application in automobile 
windshields (laminated glass), bimetal thermostats, and plywood. In these cases homogenous isotropic layers of materials are bonded 
together to form non-homogenous composite laminates. 

Composite material can be defined as a multi phase material from a combination of materials, differing in form, which remain 
bonded together, but retain their identities and properties. The composites don’t merge completely. They maintain an interface 
between each other and act in concert to provide improved properties not obtainable by any of the original components acting singly. 
The properties that can be improved by forming a composite material include strength, stiffness, corrosion resistance, wear resistance, 
attractiveness, reduction in weight, fatigue life, temperature dependant behavior, thermal insulation thermal conductivity, acoustical 
insulation etc. More recently fibre reinforced resin components that have high strength to weight and stiffness to weight ratios have 
become important in weight sensitive applications such as aircraft and space vehicles. 
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7.2. Laminated composites 

Consists of layers of at least two different materials that are bonded together. Lamination is used to combine the best aspects 
of the constituent layers in order to achieve a more useful material. The properties that can be emphasized by lamination are strength, 
stiffness, low weight, corrosion resistance, wear resistance, thermal insulation etc. Bimetals, clad metals, laminated glass, plastic based 
laminates and laminated fibrous composites fall in this category. Laminated fibrous composites are a hybrid class of composites 
involving both fibrous composites and lamination techniques. A more common name is laminated fibre reinforced composites. Here, 
layers of fibre reinforced material are built up with the fibre directions of each layer typically oriented in different directions to give 
different strengths and stiffness in the various directions. 

7.2.1. Lamina and laminate: 

A lamina is a flat arrangement of unidirectional fibres in a matrix. The fibres and the matrix do not normally dissolve 
cohesively together but do contribute to a synergetic property change. The load bearing characteristics, the strength, and stiffness are 
virtually provided by these reinforced fibres while the matrix forms the body and grips or holds the reinforcement of the fibre together. 
The function of the matrix is to support and protect the fibres and to provide a means of distributing load among and transmitting load 
between the fibres. And also it provides resistance to crack propagation and damage and determines the overall service temperature 
limitations of the composite. Fibres generally exhibit linear elastic behaviour and resinous matrix materials are visco-elastic if not 
visco-plastic. Fibre reinforced composites such as boron epoxy and graphite epoxy are usually treated as linear elastic materials since 
the fibres provide the majority of the strength and stiffness. Refinement of that approximation requires consideration of some form of 
plasticity, viscoelasticity, or both. Very little work has been done to implement those idealizations of composite material behaviour in 
structural applications. Hence in present study, fibre reinforced composites have been treated as linear elastic materials. 
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A laminate is a stack of laminae with various orientations of principal material directions in the laminae. The layers of a 
laminate are usually bound together by the same matrix material that is used in the laminae. Laminates can be composed of different 
plates of different materials or, layers of fibre-reinforced laminae. The major purpose of lamination is to tailor the directional 
dependence of strength and stiffness of a material to match the loading environment of the structural element. 
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7.3. Mathematical Modelling of Composite I-Beam 


7.3.1. Kinematics 

The theoretical developments presented in this thesis require three sets of coordinate systems, which are mutually interrelated. 
The first coordinate system is the orthogonal Cartesian coordinate system (x, y, z), for which the x and y axes lie in the plane of the 
cross-section and the z-axis parallel to the longitudinal axis of the beam. The second coordinate system is the local plate coordinates 
(n, s, and z) as shown in Fig. 7.1, wherein the n-axis is normal to the middle surface of a plate element, the 5-axis is tangent to the 
middle surface and is directed along the contour line of the cross-section. The (n, s, z) and (x, y, z) coordinate systems are related 
through an angle of orientation as defined in Fig. 7.1. The third coordinate set is the contour coordinate s along the profile of the 
section with its origin at any point O on the profile section. Point P is called the pole axis, through which the axis parallel to the z- 
axis is called the pole axis. 

To derive the analytical model for a thin-walled composite beam, the following assumptions are made: 

1. The contour of the thin wall does not deform in its own plane. 

2. The shear strain y of the middle surface is zero in each element. 

3. The Kirchoff-Law assumption in classical plate theory remains valid for laminated composite thin-walled beams. 


5 




Fig. 7.1. Definition of coordinates in thin-walled open section 


According to assumption 1, the mid surface displacement components u , v at a point A in the contour coordinate system can 
be expressed in terms of a displacements U, V of the pole P in the x, y directions, respectively, and the rotation angle (j) about the pole 


axis 


u (s,z) =U(z) sin 0(s) -V(z) cos 0(s) - ^(z)q(s), 


(7.1a) 


v(s,z) = U(z) cos 0(s) +V(z) sin P(s) + ^(z)q(s). 


(7.1b) 


These equations apply to the whole contour. The out-of-plane shell displacement w can now be found from assumption 2 (also 
known as Vlasov assumption). For each element of middle surface 
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(7.2) 


7 


dv dw 

— + — =0 


dz ds 

After substituting forv from Equations (7.1a) and (7.1b), and considering the following geometric relations: 
dx =ds cos0, (7.3) 

dy =ds sin0, (7.4) 

Eq. (7.2) can be integrated with respect to s from the origin to an arbitrary point on the contour 
w (s, z) =W(z) - U 1 (z)x(s) - V 1 (z)y(s)-^(z)w(s) , (7.5) 

Where differentiation with respect to the axial coordinate z is denoted by primes ('); W represents the average axial displacement of the 
beam in the z-direction; x and y are the coordinates of the contour in the (x, y, z) coordinate system; and 0) is the so-called sectorial 
coordinate or warping function given by 

a>(s) = J r(s)ds (7.6) 


The displacement components u, v, w representing the deformation of any generic point on the profile section are given with respect 
to the mid surface displacements u , v , w by assumption 3. 

u(s, z, n) = u (s, z), (7.7a) 


v (s, z, n) = v (s, z) - n 


du(s,z) 


ds 


(7.7b) 


, \ / \ du(s,z ) 

w (s, z, n) = W (s,z) -n—--, 

dz 

The strains associated with the small-displacement theory of elasticity are given by 


(7.7c) 
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=£ s +nk ,, 


(7.8a) 


£- = £ - + nk 

z z z ’ 


nk 


Where 


dv 


dw 


ds 


dz 


(7.8b) 

(7.8c) 


(7.9a) 


d~u - 


d~u 


d 2 u 


(7.9b) 


ds~ " dz ds dz 

All the other strains are identically zero. In equations (7.9a) and (7.9b), £ z ,k y are assumed to be zero, and s k and k sz are mid 

surface axial strain and biaxial curvatures of the shell, respectively. The above shell strains can be converted to beam strain 
components by substituting Equations. (7.9a) and (7.9b) 


£ = £° z +xk y +yk T +tf>k (0 , (7.10a) 

k z = k y sin 0 - k v cos 6 - q, (7.10b) 
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(7.10c) 



Where e z , k ^ , k y , k io and k v , are axial strain, biaxial curvatures in the x and y 

directions, warping curvature with respect to the shear center, and twisting curvature in the beam, respectively, defined as 


= W’, 

(7.11a) 

5 : 

> 

1 

II 

* 

M 

(7.11b) 

II 

l 

(7.11c) 

e 

II 

1 

■s. 

(7.lid) 

k sz =2^'. 

(7.lie) 


The resulting strains can be obtained from Equations. (7.8a)- (7.8c) and (7.10a)- (7.10c) as 

£, = £? + (x + n sin 0)k y +(y-n cos 0) k x + (&>- nq) k w , (7.12a) 

r«=nk„. (7.12b) 

7.3.2. Variational formulation 
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The total potential energy of the system can be stated as 

n= ~\\( (7 z £ z +(T zsr Z s)dv . (7.13) 

V 

the variation of the total potential energy is calculated by substituting Equations (7.12a) and 
(7.12b) into equations (7.13) 

*n= JV . [8s 0 , + (x + n sin 0)8k, +(y-n cos 0)Sk x + (o) - nq)Sk m )] + <J sz n§k sz }dv 

= +M.Sk.+M,SkJdz (7.14) 

Where N _. M c , M and M (0 are axial force, bending moments in the x and y directions, and warping moment (bi moment) with 
respect to the centroid, respectively, defined by integrating over the cross-sectional area A as: 


N, = [ a.dsdn , 

z J A z 

(7.15a) 

M f = £ <t, (y - n cos 0)dsdn , 

(7.15b) 

M = [ adx + nsin 0)dsdn, 

y JA 4 

(7.15c) 
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M m = f aAco - nq)dsdn, 

JA 

(7.15d) 


M ,= [ crndsdn, 

Ja zs 

(7.15e) 


The kinetic energy of the system is given 



T = ^- |/7(m 2 +V 2 +W 2 )dv, 

(7.16) 


Where p is a density. 

The variation of the kinetic energy is expressed by substituting the assumed displacement 

field into Eq. (7.16) as 

St = jp{USU + VSV + WSW + (q 2 + r 2 - 2rn + n 2 ) (f)S<j) 


+ ((pSU + US<j) ) [n cos 0 - (y-y p )] 



+ ( tf>SV + V8(j) ) [n cos 0 - (x-x p )] } dv 

(7.17) 


in Equation(7.14), the following geometric relations are used (Fig. 7.1): 

x - x = q cos 6 + r sin 0 

(7.18a) 
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y - y = q sin 0 - r cos 6 

In order to derive the equations of motion, Hamilton’s principle is used 


(7.18b) 


^f 2 (r-n)dt = 0. (7.19) 

substituting Equations into Equations (7.19), the following weak statement is obtained: 

0 = [" f {m 0 WSW+ [m 0 U + (m c +m 0 y )<f>]SU 
Jo ‘ 

+ [ m 0 V + (m s -m o x p )0]SV 

+ [ (in + m 0 y p )U + (m s - m 0 x p )V l 

+ ( m p +m 2 - 2m /o )(/)\S(/) 

-N Z SW' + M y SU" + M x 5V" + M w S</>" 

-2M t S</>'} dzdt. J (7.20) 

In Equation (7.20), m 0 , m c , m s , m n) , m p , m 2 are inertia coefficients, respectively, defined by 
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m 0 = I 0 j\&, 

(7.21a) 

m f = cjcos# ds, 

(7.21b) 

m s = I , J sin 6? ds. 

(7.21c) 

= 1 1 j/ ds, 

(7.21d) 

m p = !o{(« 2 +r2 ) ds ’ 

(7.21e) 

m 2 = I 2 j" ds, 

(7.2 If) 

where 


(1 0 ,1!, 1 2 ) = f p(l,n,n 2 )dn. 

Jn 

(7.22) 


7.3.3. Constitutive equations 

The constitutive equations of a k th orthotropic lamina in the laminate coordinate system are given by 
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(7.23) 


Qu Q\6 

016 066 

where Q jj are transformed reduced stiffness[12]. Axial force N _ can now be 
expressed with respect to the generalized strains (s°, k y , k , k (0 , k sz ) by 
combining Equations. (7.15a), (7.23) and (7.12a), (7.12b) 

N - = £ {q, j [s z +(x + n sin 0)k y +(y~n cos 0)k x +{oj- nq)k 0) ] + Q l6 nk sz } . (7 . 24) 

Similarly, the other stress resultants (M, M, M, M) can also be written in terms of the 

generalized strains. Consequently, the constitutive equations for a thin-walled laminated composite are obtained as 




X' 


E n E 12 E l 3 E l 4 E l5 


k°l 

M y 


^22 E 24 E 25 


K 

M x 

. = 

E 23 E m E 35 

< 

K 

M co 


E 44 E 45 


k co 

M, 


sym, E 55 


k s Z . 


Where E t] are stiffness of the thin-walled composite, and can be defined by 


(7.25) 
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E n = 

A, | ds, 

Js 

(7.26a) 

E n 

| (Aj jX + B { j sin 0)ds, 

(7.26b) 

I - 

^ 13 

|(Aj l y-B n cos 0)ds, 

(7.26c) 

E \ 4 

J(- A n 6)-B n q}ls, 

(7.26d) 

I = 

j (^16 )*, 

(7.26e) 

F - 

‘-'22 

| (Aj,x 2 + 2 B { [X sin 6 + D { 1 sin 2 d)fis , 

(7.26f) 

F - 

‘-‘22 

| (A, X xy + j (y sin 6 1 - xcos 6)-D n sin 6>cos 

(7.26g) 

F - 

‘-'24 

| (A, jX&> + 2? n (a sin 6 -qx ) - sin 6>)ds, 

(7.26h) 

F - 

‘-'25 

| (B l6 x + D l6 sin 0)ds, 

(7.26i) 

F - 

‘-'22 

J (a u j 2 -2B u ycos0 + D u cos 2 0)ds, 

(7.26j) 

F - 

‘-'24 

J (A u ya>-B n (co cos 6 + qy) + D n q cos &)ds. 

(7.26k) 

F - 

‘-'25 

\ (B l6 y-D 16 cos 0)ds, 

(7.261) 
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E A 4 = \(a u O) 2 - 25, l coq + D { { q 2 }ls. 

(7.26m) 

E 45 — j (B i6 ct) — D [6 q)ds, 

(7.26n) 

E 55 — jD 66 ds. 

(7.26o) 

A .., B .., and D .. matrices are extensional, coupling and bending stiffness, respectively, 

defined by 

( A ij , B y, D tj ) = \Q ij (hn,n 2 )dn 

(7.27) 

It appears that the laminae stiffness E jj depend on the cross-section of the composite, 

and the explicit expressions for I- section are given as follows: 

11 

3* 

(7.28a) 

^12 ~ ^11^3’ 

(7.28b) 

En = A“b a y a -B“b a , 

(7.28c) 

d 

II 

(7.28d) 
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- Kh- 


(7.28e) 

£ 2! = A“§+o, j A. 


(7.281) 

^23 O, 


(7.28g) 

£ 2 4 = k*-A“,]§. 


(7.28h) 

^25 = ^ 16 ^ 3 ’ 


(7.28i) 

£ u =te-2Br*+B,",’ 

h + n A ' v 

(7-28j) 

£ = i D ’ 

34 i2 


(7.28k) 

£,5 = s“2A-orA. 


(7.281) 

£4 4 = te-2B,>„ + A"' 

12 12 11 

(7.28m) 

£ 45 = 0 , 


(7.28n) 

E 55 = D^d k , 


(7.28o) 
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y 



Fig.7.2. Geometry of a thin-walled doubly symmetric 
Composite I-beam. (bi=b 2 and ti=t 2 ) 

In the above equations, repeated index denotes summation. Index k varies from 1 to 3 
whereas a varies from 1 to 2 implying top and bottom flanges (1,2) and web (3) as 
shown in fig. b k denotes width of flanges and web. 
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Similarly , the explicit expressions of inertia coefficients for I-section are given by 


m o = 7 X> 
>n c =I?b a , 


(7.29a) 

(7.29b) 


m s = I fa. 


m„ = K 


— + y 2 & 

* a a 


+ I > b A 

0 12 


m 2 =l 2 b k . 


(7.29c) 

(7.29d) 

(7.29e) 

(7.29f) 


7.3.4. Equation of motion 

The equation of motion of the present study can be derived by integrating the derivatives of the varied quantities by parts and 
collecting the coefficients of Su, Sv, Sw,andSt/> . 

N' z =m 0 W, (7.30a) 

M" y = m o U + (m c + m D y p )l>, (7.30b) 
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M x = m 0 V + (m s + m o x p )o, (7.30c) 

M a + 2M, t =(m c +m 0 y p )U+(m s -m 0 x p )V+{m p +m 2 -2m a )0. (7.30d) 

The natural boundary conditions are of the form: 

SW: N z , (7.31a) 

SU:M' y , (7.31b) 

SU : M y , (7.31c) 

SV:M X , (7.3 Id) 

SV :M x , (7.3le) 

d®:M C0 +2M t , (7.3 If) 

&5>:M (0 , (7.31g) 
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The above equations are most general form for flexural-torsional vibration of a 
thin-walled laminated composite beam with an I-shaped section, and the dependent variables, u, v, w and (j> are fully coupled. If the 

stacking sequence of the web is symmetric, E u = E u = 0, and the thin-walled composite is symmetric with respect to z-axis, 
E n = E 15 = E 24 = 0. Further, if both the web and flange are balanced laminates ( ± 0 pairs of layers) A. b = = 0, and thus, 

E 25 = 0. Finally, Eq. (7.30a) can be simplified to the uncoupled differential equations as: 


( EA) com W"=pAW , 

(7.32a) 

~(EI y ) com U w = pAU, 

(7.32b) 

-( EI x ) com V' v = pAV , 

(7.32c) 

-(EIJ com ® iV +(GJ) eom = plA 

(7.32d) 

From the above equations, ( EA) com represents axial rigidity, (El x ) com 

and (EI y ) com 

represent flexural rigidities with respect to x and y axes, ( EI (0 ) com and 

walled composite, respectively, written as 

(GJ) com represent warping and torsional rigidities of the thin 

(EA)com = Ell, 

(7.33a) 

(Ely)com = E22, 

(7.33b) 

(Elx)com = E33, 

(7.33c) 
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(EIoo) com — E 44 , 


(7.33d) 


(GJ)com =4 E 55 , (7.33e) 

In Eqs. (7.32a)-(7.32d), Idenotes the polar moment of inertia. It is well known that the three distinct vibration modes, flexural 

vibration in the x and y directions, and torsional vibration, are identified in this case, and the corresponding vibration frequencies are 
given by the closed form for simply supported boundary conditions[3]: 



(7.34a) 

(7.34b) 

(7.34c) 


(7.34d) 

where a> z , a> x , co y , a> 0 are axial frequency, flexural frequencies in the x and y 
directions, the torsional vibration frequency, respectively 
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7.3.5. Finite Element Model 


The present theory for thin-walled composite beams described in the previous section was implemented via a displacement based 
finite element method. The generalized displacements are expressed over each element as a linear combination of the one¬ 
dimensional Lagrange interpolation function <p ; and Hermite-cubic interpolation function <p t associated with node j and the nodal 
values: 


n 


w =2 >#,■’ 
j= i 

(7.35a) 

n 

U^UjVj, 

7=1 

(7.35b) 

II 

iM 5 

< 

(7.35c) 

7=1 

(7.35d) 


23 



Substituting these expressions into the weak statement in Eq. (7.20), the finite element 
model of a typical element can be expressed as a standard Eigen value problem 
([AT]-A[Af]{A} = {0}, (7.36) 


where [K] is the element stiffness matrix 


*11 

K ,2 

K 13 

K 14 


k 22 

K 2i 

k 2 , 



K 

*34 

sym . 



^44 


and [M] is the element mass matrix 


M ii 

M 12 

M 13 

M u 


^22 

^23 

Mu 



M 33 

Mu 

sym . 



M 44 


(7.37) 


(7.38) 


The explicit form of [K] and [M] are given by 

K" =^E n .i//\i//'j dz, (7.39a) 

K-j 2 = £ E xr y/\ i/'jdz, (7.39b) 
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</ = ££ l3 .^>'' ;t fe, (7.39 c) 

K^=\ l 2E l5 yyjdz, (7.39d) 

^ 2 = ['£ 22 y>".*, (7.39e) 

^ 24 = ££ 24 y>".*, (7.390 

^ 3 = ££ 33 y>'V^, (7.39g) 

= £(^y, ¥ "-2E, i¥ '\ ¥ ' } )dz, (7.39h) 

^ 4 = £(£ 44 ^", ^" 2 ^ 55 .K- VSj )dz, (7.39i) 

M,y = G,f = M 33 = £m 0 ^ jdz, (7.40a) 

= £ o (w c + m 0 y p )i// l i// J dz, (7.40b) 

M 34 = £ (m s - m n x p ) ¥¥j dz, (7.40c) 
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all the other components are zero. 

In equations (7.36), { A } is the eigenvector of nodal displacements corresponding to an 
Eigen value {a} = {wuv</)} t 


(7.41) 
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Chapter8 


EFFECT OF VARIOUS PARAMETERS ON VIBRATIONS OF UNIFORM 

COMPOSITE 

I-BEAM 

A thin-walled composite beam with I-section and length L = 8 m is considered in order to investigate the effects of various 
parameters namely fibre orientation, modulus ratio, height-to-thickness ratio and boundary conditions on the natural frequencies and 
mode shapes. The geometry of the I-section is shown in Fig. 8.1 and the following engineering constants are used for finite element 
analysis. 
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8.1. Material properties 


Graphite-Epoxy is used for composite beam analysis. The properties of Graphite-Epoxy are 

• Modulus of Elasticity Ej =1.38* 10 11 N/m 2 , 

• Density p = 1590 Kg/m 3 , 

• EJEi = 40, 

• G 12 /E 2 = 0.6, v 12 = 0.25 Where subscripts 1 and 2 indicate fibre direction and perpendicular to fibre direction, 
respectively. 

For convenience, the following non-dimensional natural frequency is used: 

col 2 

m= - 

b, 



8.1.1. Element used 

The element considered for the composite beam analysis is Shell99. Shell99 is used for layered applications of a structural 
shell model. The element has six degrees of freedom at each node, translations in the nodal x, y, and z directions and rotations 
about the nodal x, y, and z-axes. The element is defined by eight nodes, average or comer layer thicknesses, layer material 
direction angles, and orthotropic material properties. It usually has a smaller element formulation time. Shell99 allows up to 250 
layers. If more than 250 layers are required, a user-input constitutive matrix is available. The figure 8.2 shows Shell99, Linear 
Layered Structural Shell. 
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Figure 8.2.SHELL99, Linear Layered Structural Shell. 


8.2. Modal analysis 

A composite I- beam is modeled using ANSYS with different fibre orientation and taking 20 layers are laid one over another 
and modal analysis is carried out . The effect of fibre orientation, elastic modulii and height to thickness ratio on natural 
frequencies and mode shapes of composite I-beam is studied. 
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8.2.1 Effect of Fibre angle rotation in top and bottom flanges 


First, the top and bottom flanges are considered as angle-ply laminates [9/ -9], and the web laminate is assumed to be 
unidirectional in I-beam. In addition, the stacking sequence of the beam is symmetric with respect to x-axis. The coupling 
stiffnesses E 12 , E 13 , E 15 , E 24 and E 25 become zero, but E 35 does not vanish due to unsymmetric stacking sequence of each flange. 

Accordingly, flexural vibration in the y-direction is uncoupled, whereas the flexural vibration in the x-direction and the 
torsional vibration are coupled. The lowest four non-dimensional natural frequencies by the finite element analysis and the 
orthotropic closed-form solutions, which neglects the coupling effects obtained from equations (7.34a) to(7.34d) for each mode 
are given in Table 8.1.For unidirectional fibre orientation, the lowest four natural frequencies by the finite element analysis exactly 
correspond to the first flexural mode in the x-direction, torsional mode, second flexural mode in the x-direction, and flexural mode in 
the y-direction by the orthotropic closed-form solution, respectively. As the fibre angle changes, the variation in natural 
frequencies is illustrated in Fig. 8.3, and the results by finite element analysis and closed form solution show slight discrepancy. 
This is because of the coupling stiffnesses which are neglected in the closed-form analysis. 
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Table 8.1: Non-dimensional natural frequencies with respect to the 
fibre angle change in the flanges. 


^"Natural freq 
Fibre angle . 


Closed form 

FEM 

0) ; ®,2 

co x a> 2 co 3 co 4 


15.859 

19.481 

39.419 

41.085 

15.859 

19.48 

39.419 

41.085 

12.163 

25.562 

36.564 

38.504 

12.214 

22.365 

36.294 

38.042 

6.4308 

30.126 

23.386 

32.381 

6.4315 

23.146 

24.142 

31.452 

3.1165 

34.536 

16.498 

29.365 

3.1165 

13.599 

25.465 

24.143 


















8.3.Percentage variation of frequencies closed form Vs FEM Solutions 


The difference and variation between FEM and closed form solution is presented from table 8.2 to table 8.5, for a simply 
supported beam with respect to the fibre angle change in flanges. 


Table 8.2. Percentage variation in first natural frequencies of a simply supported 
composite beam with respect to fibre angle change in flanges 


Fibre 

Closed form 

FEM 

difference 

% 

Angl 

Natural 

Natural 


varia 

e 

Frequencies 

Frequencies 


tion 



a> l 



0 

15.859 

15.859 

0.0 

0.0 

15 

12.163 

12.214 

-0.051 

0.42 

30 

6.4308 

6.4315 

-0.0004 

0.01 

45 

3.1165 

3.1165 

0.0 

0.0 
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J 










Oh 


0 

19.481 

19.48 

0.001 

15 

25.562 

22.365 

3.197 

30 

23.126 

23.146 

-0.02 

45 

34.536 

33.599 

0.937 

60 

30.052 

25.946 

4.106 

45 

20.542 

20.946 

-0.404 

90 

14.512 

14.512 

0.0 


0 


12.5 

0.086 

2.713 

13.66 

1.966 


0 





Table 8.4. Percentage variation in third natural frequencies of a simply supported 
composite beam with respect to fibre angle change in flanges 


Fibre 

Closed form 

FEM 

Difference 

%Variation 

Angle 

Natural 

Natural 




frequencies 

Frequencies 





<X> 3 



0 

39.419 

39.419 

0.0 

0.0 

15 

36.564 

36.294 

0.27 

0.738 

30 

23.386 

24.142 

-0.756 

3.232 

45 

26.498 

25.465 

1.033 

3.898 

60 

14.635 

15.415 

-0.78 

5.329 





45 


14.546 


19.512 


-4.966 


0.341 


90 


10.562 


10.562 


0.0 


0.0 





Table 8.5. Percentage variation in fourth natural frequencies of a simply supported 
composite beam with respect to fibre angle change in flanges 


Fibre 

Closed form 

FEM 

Difference 

% 

Angle 

Natural 

Natural 


Variation 


Frequencies 

Frequencies 




(O y 

0) 4 



0 

41.085 

41.085 

0.0 

0.0 

15 

38.504 

38.042 

0.462 

0.011 

30 

32.381 

31.452 

0.929 

0.0286 

45 

29.365 

24.143 

2.222 

0.045 

60 

24.912 

24.102 

0.80 

0.0286 





45 


26.546 


26.451 


0.095 


0.0035 


90 


26.469 


22.643 


3.826 


0.1445 






fiber angle ,in degrees 
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Fig.8.3 .Variation of the non-dimensional natural frequencies of a simply 
supported composite beam with respect to fibre angle change in 
flanges 

The variation of natural frequencies of a simply supported composite beam by orthotropic closed-form solution with fibre 
angle change in the flanges for the lowest four vibration mode shapes are shown in Fig.8.3. It is seen that the mode shape 
corresponding to the lowest natural frequency is the first flexural mode in z-direction ( co x] ), which is well below the other three types 

of modes. In general, co xl , co x2 and co y decrease as the fibre angle increases. In torsional vibration, however, the laminate torsional 

rigidity, {GJ} CO m, has an important effect on vibration characteristics. Thus, placing the fibre angle at ±45° leads to considerable 
improvement in the torsional vibration frequency as can be seen for the compressive load-carrying capacity in the in the paper [11]. 

8.4. Effect of Fibre angle rotation in Web: 

The fibre angle is rotated in the web instead of the flanges for the same composite I-beam referred in fig.8.1. The 
variations of natural frequencies for various modes by finite element solution with respect to the fibre angle variation in the 
web are shown in Fig. 8.4. 

Table 8.6: Nondimensional natural frequencies with respect to the fibre 
angle change in the web. 


43 




Closed form 

FEM 

Natural 









frccV 






co 2 

<X>3 

0) 4 

Fibre \ 









angle \ 









0 

15.859 

19.481 

39.419 

41.085 

15.859 

19.481 

39.419 

41.085 

15 

15.494 

21.68 

39.614 

38.421 

15.804 

20.968 

38.83 

39.524 

30 

15.483 

26.296 

39.564 

36.689 

15.491 

24.148 

36.542 

38.562 

45 

15.446 

24.953 

39.463 

35.569 

14.158 

25.053 

35.231 

34.024 

60 

15.448 

24.804 

39.452 

35.365 

14.058 

24.481 

30.941 

35.21 

45 

15.485 

21.126 

39.253 

35.286 

14.022 

22.142 

29.146 

34.462 

90 

15.499 

19.458 

38.984 

35.252 

13.966 

14.165 

28.125 

33.912 


Table 8.7. Percentage variation in first natural frequencies of a simply supported 
composite beam with respect to fibre angle change in web 
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Fibre 

Angle 

Closed form 

Natural Frequencies 

FEM 

Natural 

Frequencies 

Difference 

% Variation 


®vl 

co x 



0 

15.859 

15.859 

0.0 

0.0 

15 

15.494 

15.804 

0.31 

0.02 

30 

15.483 

15.491 

0.008 

0.0005 

45 

15.446 

14.158 

1.288 

0.083 

60 

15.448 

14.058 

1.39 

0.0589 

45 

15.485 

14.022 

1.463 

0.0944 






90 


15.499 


13.966 


1.533 


0.0989 


Table 8.8. Percentage variation in natural frequencies of a simply supported 
composite beam with respect to fibre angle change in web 



Closed form 

FEM 



Fibre 

Angle 

Natural 

Frequencies 

Natural 

Frequencies 

Difference 

% 

Variati 

a o 

co 2 


on 


0 

19.481 

19.481 

0.0 

0.0 

15 

21.68 

20.968 

0.712 

0.032841 

30 

26.296 

24.148 

2.148 

0.081685 


6 








45 


24.953 


25.053 


60 24.804 24.481 

45 21.126 22.142 

90 19.458 14.165 


0.1 


0.004008 


0.323 

1.016 

5.293 


0.013022 

0.048092 

0.272022 


4 


■7 






Table 8.9. Percentage variation in natural frequencies of a simply supported 
composite beam with respect to fibre angle change in web 



Closed form 

FEM 



Fibre 

Angle 

Natural 

Frequencies 

Natural 

Frequencies 

Difference 

% 

Variation 



<» 3 



0 

39.419 

39.419 

0.0 

0.0 

15 

39.614 

38.83 

0.784 

0.019791 

30 

39.564 

36.542 

3.022 

0.076383 

45 

39.463 

35.231 

4.232 

0.10724 





60 


39.452 


30.941 


8.511 


0.215731 


45 

39.253 

29.146 

10.107 

0.257484 

Table 8.10. Percentage variation in 

natural frequencies of a simply 

90 

38.984 

28.125 

10.859 

0.27855 

supported composite beam with 

respect to fibre angle change in web 



Closed form 

FEM 



Fibre 

Angle 

Natural 

Frequencies 

Natural 

Frequencies 

Difference 

% 

Variation 



0) 4 



0 

41.085 

41.085 

0.0 

0.0 
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15 


38.421 


30 

45 

60 

45 


36.689 

35.569 

35.365 

35.286 


90 


35.252 


39.524 


-1.103 


0.028 


38.562 0.124 0.003 

34.024 -1.455 0.039 

35.21 0.155 0.004 

34.462 -0.0146 0.004 

33.912 1.34 0.038 


1 




45 

40 
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15 


30 


45 


60 


75 


90 


fiber angle, in degrees 
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Figure.8.4.Variation of the non-dimensional natural frequencies of a 

simply supported composite beam with respect to fibre angle 
change in the web 

It is seen that the natural frequencies corresponding to the first and second flexural modes in the z-direction ( co xl , co x2 ) are 

almost invariant for entire range of fibre angle variation in the web. The torsional frequency reaches its maximum near 0 = 45° 
showing similar trends as for the previous case. 

In order to investigate the coupling effect further, the top flange is considered as an angle-ply laminate [ 0/ — 0], while the 
bottom flange and web are unidirectional. In this case, the thin-walled composite beam is unsymmetric with respect to z-axis. As the 
fibre angle changes, the results by the finite element analysis and the orthotropic closed form solution show significant discrepancy 
implying that the coupling effects in this case become no longer negligible. 

8.5. Mode shapes 

The mode shapes corresponding to the first four lowest modes with unidirectional fibre and with fibre angle 30° in the top 
flange are illustrated in Figs. 8.5 and 8.6.It can be readily seen that the beam with fibre angle 30° in the top flange exhibits strong 
flexural-torsional coupling. In fact, the coupling stiffnesses En and E 25 vanish while E n ,E l5 ,E 24 and E35 do not vanish due to globally 
unsymmetric lay up. That is, the orthotropic solution is no longer valid for unsymmetric laminated beams, and fully coupled 
flexural-torsional vibration should be considered for accurate analysis. 
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Fig.8.5. Mode shapes of the composite beams with unidirectional fibre. 
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(a) Mode 1 (b) Mode 2 









( c) Mode 3 (d) Mode 4 



Fig.8.6. Mode shapes of the composite beams with fibre angle 30° in the web. 
















8.5.1. Effect of Elastic modulus ratio 

The next case study shows the effects of modulus ratio (E l /E 2 ) of composite beams on the natural frequencies and mode 
shapes for simply supported and clamped composite beams (Figs. 8.7 and 8.8). The stacking sequence of the top and bottom flanges 
are [0/90] s , and web is unidirectional. It is observed that the natural frequencies increase with increasing orthotropy ( E l /Ei) for 
both simply supported and fixed boundary conditions. 

8.6. Effect of Height-to-thickness ratio: 

The effects of the height-to-thickness ratio (b 3 /t) of thin-walled composite beams on the natural frequencies and mode shapes 

are illustrated in Fig. 8.7 and fig. 8.8 for simply supported and clamped boundary conditions. Torsional frequency is found to be very 
sensitive to height-to-thickness ratio while the other frequencies almost invariant for both boundary conditions. 
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Fig.8.7.Variation of the non-dimensional natural frequencies of a simply Supported composite beam with respect to modulus ratio 
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Fig.8.8.Variation of the non-dimensional natural frequencies of fixed 
composite beam with respect to modulus ratio 
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8.6.1. Effect of Height-to-thickness ratio 

The variation of the non-dimensional natural frequencies of a simply supported composite beam and fixed beam with respect to height 
-to-thickness ratio are shown in figure 8.9 and 8.10.. 
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Fig.8.9.Variation of the non-dimensional natural frequencies of a simply supported 
composite beam with height-to-thickness ratio. 
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Fig.8.10.Variation of the non-dimensional natural frequencies of a fixed composite 

beam with height-to-thickness ratio. 
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8.6.2. Harmonic analysis. 

The harmonic analysis is carried out for a cantilever doubly symmetric composite I-beam by applying a load of 1 Newton at free 
end. The analysis is carried out by varying fibre angle rotation in top and bottom flanges of I- beam, elastic modular ratio and height to 
thickness ratio of I- beam and results are obtained and presented in the tabular form. These results are also depicted in the graphical 
form as Frequency Vs. Amplitude in figures 8.11, 8.12 and 8.13. 

8.6.2.1. Effect of fibre angle rotation in top and bottom flanges 

In this example the top and bottom flanges are considered as angle-ply laminates [9/ -9], and the web laminate is assumed 
to be unidirectional. The numbers (0-90) shown in graph are fibre angle rotations [9/ -9] in degrees. 
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Amplitude,in metres 



Fig.8.11. Frequency Vs. Amplitude for fibre angle rotation in two flanges. 
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From the Fig.8.11, for unidirectional fibre (zero fibre angle rotation), the amplitude of vibration is maximum at 10 Hz and 
decreases suddenly and slowly increases after wards. For a 45 degree fibre angle rotation, the value of the amplitude increases from 
20-30 Hz and decreases afterwards. 

S.6.2.2. Effect of Elastic modulus ratio 

This analysis shows the effects of modulus ratio (E l /E 2 ) of composite beams on the harmonic response for cantilever 
composite I-beams .The stacking sequence of the top and bottom flanges are [0/90] s , and web is unidirectional is considered in this 
analysis The figure 8.12 shows the variation of natural frequencies Vs amplitude for elastic modular ratios (E 1 /E 2 ) from lto 40. 
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Fig.8.12. Frequency Vs. Amplitude for various elastic modulus ratios 

From the fig.8.12, it is clear that as the frequency increases, the value of the amplitude of vibration for a cantilever composite 
beam follows sinusoidal curve, for all modular ratios from 10-40 , except for modular ratio one and the values of amplitude of 
vibration decreases in frequency step 10-20 Hz and gradually increases for next step 20-30Hz and then decreases. 
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8.6.2.3. Effect of Height-to-thickness ratio 

The effects of the height-to-thickness ratio (b 3 /t) of thin-walled composite beams on the Harmonic response illustrated in 
Fig.8.13 for cantilever boundary conditions for the different values of height-to-thickness ratio (b 3 /t) of 10 to 40. 
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Fig. 8.13. Frequency Vs. Amplitude for various beam height to thickness ratios 
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From the fig. 8.13, it shows that the amplitude of vibration is maximum for frequency of 40 Hz. for a height to thickness ratio of 10. 
For other height to thickness ratios, it is almost constant for all other values of frequencies. It clearly shows that the amplitude of 
vibration decreases by increasing the height to thickness ratio at all frequencies. 
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Chapter 9 

MODAL AND HARMONIC ANALYSIS OF TAPERED COMPOSITE I- 

BEAM WITH A SMALL END FIXED 
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Chapter 9 


MODAL AND HARMONIC ANALYSIS OF TAPERED COMPOSITE I- BEAM 

WITH A SMALL END FIXED 

A thin-walled tapered composite beam with I-section and length L= 8 m is considered in order to investigate the effects of 
fibre orientation and modulus ratio on the natural frequencies and mode shapes for a cantilever boundary condition. The Finite 
Element model of a tapered I-beam is shown in Fig. 9.1.The engineering constants and element used is same as in the case of a 
uniform I-beam, mentioned in 8.1 and 8.1.1. In this case, the thickness of the flanges and web is constant but the height and width 
of the beam is tapered to half of the actual size from one end to another. The smaller end is fixed and another end is free. For 
Harmonic analysis 1 Newton load is applied at free end of the cantilever composite I-beam and results are obtained are presented in 
graphical form. 
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Fig9.1. Tapered composite I- beam with layers 
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9.1. Modal Analysis 

The modal analysis of a cantilever I-beam is done using ANSYS and the steps involved to do modelling and analysis to get the 
natural frequencies and mode shapes is shown below in the form of flow chart. The different mode extraction methods in ANSYS are 

• Sub space 

• Block Lanczos 

• Power dynamics 

• Reduced 

• unsymmetric 

Flow chart to show the Procedure to do Modal analysis using ANSYS 

Choose the Element type 

k 

Choose the Appropriate Real Constants and Material Constants 

k 

Create key points at the start and end of the beam 

U 

Create a line between the above points 


& 

k 


Mesh the beam with the size control 
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Apply Boundary Conditions 


U 

Change the analysis type to MODAL 

Solution > Analysis Type >New Analysis >Modal 

1 

Set the Analysis Options 

Solution > Analysis options >Extract 10 mode>OK 

Enter <50> for the ending frequency 

U 

Solve the load Set 

U 

ANSYS generates a Substep Result for each Natural Frequency & Mode shape 

List RESULTS summary 

General Post Processor >List Results>Results Summary 
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Read Results for a Substep 


General Post Processor > Read results > First Set 

PLOT deformed geometry 
Read Results for Next Substep 

General Post Processor > Read results > Next Set 


9.1.1. Effect of fibre angle rotation in top and bottom Flanges: - 

The top and bottom flanges are considered as angle-ply laminates [6/ -6], and the web laminate is assumed to be 
unidirectional. In this case, the lowest four natural frequencies by the finite element analysis exactly correspond to the first flexural 
mode in the x-direction, and flexural mode in the y-direction, second flexural mode in the x-direction, torsional mode by the 
orthotropic closed-form solution, respectively. It is shown that from fig.9.2, the natural frequencies are increasing with increasing 
fibre angle rotation. Up to 45° fibre angle, very small variation in the natural frequencies for first three modes. 
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Fig 9.2. Variation of the non-dimensional natural frequencies of a tapered 
cantilever composite beam with small end fixed with respect to fibre angle change in flanges 
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9.1.2. Mode shapes 

The mode shapes corresponding to the first four lowest modes with unidirectional fibre in web and fibre angle 45° in the top 
and bottom flanges are illustrated in Figs.9.3. 
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Fig. 9.3. Mode shapes of the tapered composite beam with fibre angle 45 0 in top and bottom flanges 

9.1.3. Effect of fibre angle rotation in web 


In this case the natural frequencies for first mode are constant for all fibre angles as shown in fig.9.4. For second and third modes 
frequency is constant from fibre angle 0- 60° and after that a small increment is observed. For fourth mode, natural frequency is increased 
up to fibre angle 60° and then decreases. 



















Fig. 9.4. Variation of the non-dimensional natural frequencies of a Tapered Cantilever 
composite beam with small end fixed with respect to fibre angle change in web. 
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9.1.4. Effect of modular ratio 



Fig.9.5. Variation of the non-dimensional natural frequencies of a Tapered Cantilever composite beam with small end fixed with 

respect to fibre angle. 
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The next analysis shows the effects of modulus ratio (EJ E 2 ) of composite beams on the natural frequencies and mode 
shapes for tapered Cantilever composite I- beam. The stacking sequence of the top and bottom flanges are [0/90] s , and web is 
unidirectional is considered for analysis. It is observed that the natural frequencies increase with increasing orthotropy ( E t /E 2 ). 
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9.2. Harmonic analysis 

Any sustained cyclic load will produce a sustained cyclic response in a structural system. In order to predict the sustained dynamic 
behaviour of a structure, Harmonic response analysis can verify whether or not a designed structure will successfully overcome resonance, 
fatigue and other harmful effects of forced vibration. 

The harmonic response analysis was applied to a composite cantilever beam and its structural response was investigated. A load of 
1.0 N at a frequency 0-50 Hz. is applied at free end and was used to examine the composite beam to laminates stacking. The pulse was 
localised at the centre node in the free edge. The study is carried out by varying fibre angle rotation in top and bottom flanges of I- 
beam and elastic modular ratio and results are obtained and represented in the graphical format as Frequency vs. Amplitude. The 
detailed steps involved in harmonic analysis of a beam using ANSYS are shown in flow chart form. 

Flow chart to show the Procedure to do Harmonic analysis of a Cantilever beam using ANSYS 

Choose the Element type 

4 - 

Choose the Appropriate Real Constants and Material Constants 

4 

Create key points at the start and end of the beam 

4 - 

Create a line between the above points 

4 - 
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Mesh the beam with the size control 

* 

Apply Boundary Conditions 

4 

Apply 1.0 N Load at the free end 

4 

New analysis > Harmonic 

4 

Set the Analysis Options 

xx 

Set the Solution Method 

xx 

Set the DOF print out format to “Amplitude and Phase” <OK> 

♦ 

Set the Tolerance to 1.0 e -8 < OK> 

\7 

Set the Frequency Sub Steps 

\7 
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Solution > Load step Options - Time Frequency > Frequency and Sub Steps 



Set the Harmonic Frequency Range 0-50 Hz 



Set the Number of Sub steps to 100 



Set to be STEPPED 


SOLVE the model 
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9.2.1. Effect of fibre angle rotation in top and bottom Flanges 

From the fig.9.6, it is clear that as the frequency increases, the value of the amplitude of vibration for a cantilever composite beam 
follows sinusoidal curve, for all fibre angles. For fibre angle 60°, the value of amplitude of vibration is high at frequency of 5 Hz. 
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Displacement,in meters 



Fig.9.6. Frequency Vs. Displacement for fibre angle rotation in two flanges 







9.2.2 .Effect of fibre angle rotation in web 



Fig.9.7. Frequency vs. displacement for fibre angle rotation in web 








From the fig.9.7, it is seen that as the frequency increases, the value of the amplitude of vibration for a cantilever composite beam 
follows sinusoidal curve, for all fibre angles , same as above case, and for fibre angle 45°, the value of amplitude of vibration is high 
at frequency of 5 Hz . 
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Displacement, in meters 


9.2.3. Effect of modular ratio 

As shown in fig.9.8, a small variation in amplitude of vibration for all modular ratios except modular ratio of 30. For modular ratio of 
30, the minimum displacement is at frequency of 20 Hz. 
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Fig.9.8. Frequency Vs. displacement for various modular ratios 
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Chapter 10 

MODAL AND HARMONIC ANALYSIS OF TAPERED COMPOSITE I- BEAM 

WITH BIG END FIXED 
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Chapter 10 


MODAL AND HARMONIC ANALYSIS OF TAPERED COMPOSITE I- BEAM 

WITH BIG END FIXED 

In this case the same tapered beam is taken and the big end is fixed and small end is free instead of small end fixed and big end 
free as in previous case. The same engineering constants and element are used as for previous case. For harmonic analysis 1.0 N load 
is applied at free end of the cantilever beam and other end is kept free. The results obtained are presented in graphical form with 
displacement Vs natural frequency. 

10.1. Modal Analysis 

A composite cantilever I-beam is modelled with element shell99 and the material properties are taken as discussed in 8.1. Modal 
analysis is done by fixing the big end and kept smaller end free. The effect of fibre angle rotation in top and bottom flanges, elastic modulii 
and height to thickness ratio on natural frequencies and mode shapes are discussed below. 

10.1.1. Effect of fibre angle rotation in top and bottom Flanges 

The top and bottom flanges are considered as angle-ply laminates [9/ -0], and the web laminate is assumed to be unidirectional. 
In this case, the lowest four natural frequencies by the finite element analysis exactly correspond to the first flexural mode in the x- 
direction, flexural mode in the y-direction, second flexural mode in the x-direction, and torsional mode by the orthotropic closed- 
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form solution, respectively. It is clear that from fig. 10.1, the natural frequencies are increasing with increase in fibre angle 
rotation. Very small variation in the natural frequencies for first three modes is observed from 0° and up to fibre angle 45°. 
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FiglO.l. Variation of the non-dimensional natural frequencies of a tapered 
cantilever composite beam with big end fixed with respect to fibre 
angle change in flanges 


10.1.2. Mode shapes 

The mode shapes corresponding to the first four lowest modes with unidirectional fibre in web and fibre angle 45° in the top 
and bottom flanges are illustrated in fig. 10.2. 
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(a).Mode 1 



(b). Mode 2 
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(a).Mode 3 
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Fig. 10.2. Mode shapes of the tapered composite beam with fibre angle 45 0 in top and bottom flanges 
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10.1.3. Effect of fibre angle rotation in web 


In this case the natural frequencies for first mode are constant for all fibre angles as shown in fig. 10.3. For second, third and fourth 
modes frequency is constant from fibre angle 0- 45° and a small increment in each mode is observed . 



Fiber angle, in degrees 
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Figl0.3. Variation of the non-dimensional natural frequencies of a Tapered 

cantilever composite beam with big end fixed with respect to fibre angle 
change in web. 
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10.1.4. Effect of modular ratio 

The effects of modulus ratio ( E l /E 2 ) of composite beams on the natural frequencies and mode shapes for Tapered Cantilever 
composite beam is shown below. 



0 10 20 30 40 

Elastic modulus ratio 
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Fig. 10.4. Variation of the non-dimensional natural frequencies of a Taper cantilever composite beam with big end fixed with respect 

to modulus ratio. 

The stacking sequence of the top and bottom flanges are [0/90] s , and web is unidirectional. It is observed that the natural 
frequencies increase with increasing orthotropy ( E l /Ei). 
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10.1.5. Harmonic analysis 


The harmonic analysis is carried out for a tapered cantilever composite I-beam and load of 1.0 N is applied at free end. The 
study is carried out by varying fibre angle rotation in top and bottom flanges of I- beam and elastic modular ratio and results are 
obtained and represented in the graphical form as Frequency Vs. Displacement. 



Fig 10.5. Frequency vs. Displacement for fibre angle rotation in two flanges 
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Natural frequency,Hzs 


Fig. 10.6. Frequency vs. Displacement for fibre angle rotation in web 


10.1.6. Effect of fibre angle rotation in top and bottom Flanges 







From the fig. 10.5, it is clear that as the frequency increases, the value of the amplitude of vibration for a cantilever composite 
beam follows sinusoidal curve, for all fibre angles. For fibre angle 45°, the value of amplitude of vibration is high at frequency of 5 
Hz . 

10.1.7. Effect of fibre angle rotation in web 

From the fig. 10.6, it is seen that as the frequency increases, the value of the amplitude of vibration for a cantilever composite 
beam follows sinusoidal curve, for all fibre angles same as above case. 


10.1.8. Effect of modular ratio 

From the fig. 10.7, it is evident that there is a small variation in amplitude of vibration for all modular ratios except for modular 
ratio of 40. For modular ratio of 40, the minimum displacement is at frequency of 20 Hz. 
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Fig. 10.7. Frequency vs. Amplitude for various modular ratios. 
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10.2. Composite I-beam with three circular holes in web 


In this case, uniform I- beam is considered and three circular holes are made in the web at 2m apart with holes diameter of 32 
mm as shown in fig. 10.8. The beam is modeled using ANSYS with the same material properties mentioned in the chapter 8. Modal 
and harmonic analysis is done for a cantilever boundary condition. Harmonic analysis is carried out by applying a load of 1.0 N at free 
end and nodal results at free end are obtained in the form of amplitude Vs natural frequency graphs for various fibre angle rotations 
and modular ratios. The Finite element model of a uniform doubly symmetric composite I-beam with holes in web is shown in 
fig. 10.8. 
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Fig. 10.8. Uniform composite I-beam with circular holes in web 


10.2.1. Modal Analysis 

The modal analysis is carried out with the top and bottom flanges as angle-ply laminates [9/ -0], and the web laminate is 
assumed to be unidirectional in I beam. In this case, the lowest four natural frequencies by the finite element analysis exactly 
correspond to the first flexural mode in the x-direction, and flexural mode in the y-direction, second flexural mode in the x- 
direction, torsional mode by the orthotropic closed-form solution, respectively. 

10.2.1.1. Effect of fibre angle rotation in top and bottom Flanges 
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Fig. 10.9. Variation of the non-dimensional natural frequencies of a cantilever 

Composite beam with circular holes in web with fibre angle change in flanges 

From fig. 10.9, it is clear that natural frequencies for fourth mode increases with increasing fibre angle from 0 to 45° , 
maximum at 45° and then decreases afterwards. In second mode, the natural frequency at 45° fibre angle found to be minimum. 
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10.2.1.2. Mode shapes 


The mode shapes corresponding to the first four lowest modes with unidirectional fibre in web and fibre angle 45° is taken 
in the top and bottom flanges are illustrated in Fig. 10.10. 
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Fig. 10.10. Mode shapes of the Cantilever composite beam with circular holes in web with fibre angle 45 0 in top and bottom flanges 


10.2.1.3. Effect of fibre angle rotation in web 
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FiglO.l 1. Variation of the non-dimensional natural frequencies of a Cantilever composite 
beam with circular holes in web with to fibre angle change in flanges 


From fig. 10.11, it is observed that natural frequencies of third mode increases with increasing fibre angle from 0 to 45° reaches 
maximum at 45° and then decreases afterwards. For all other modes, there is no variation in natural frequencies by changing the fibre 
angle in web. 
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10.2.1.4. Effect of Elastic modular ratio 



Fig. 10.12. Variation of the non-dimensional natural frequencies of a Cantilever 
composite beam with circular holes in web with elastic modular ratio 

In this, it shows the effects of modulus ratio (E l /E 2 ) of composite beams on the natural frequencies and mode shapes for 
uniform Cantilever composite beam with circular holes in web (Fig. 10.12). The stacking sequence of the top and bottom flanges 
are [0/90] s , and web is unidirectional. It is observed that the natural frequencies increase with increasing orthotropy ( E l /Ei). 
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10.3. Harmonic Analysis 

The harmonic analysis is carried out for a cantilever composite uniform doubly symmetric I-beam and load of 1 Newton is 
applied at free end. The study is carried out by varying fibre angle rotation in top and bottom flanges of I- beam and elastic modular 
ratio and results are obtained and represented in the graphical form as Frequency vs. Displacement. 

10.3.1. Effect of fibre angle rotation in top and bottom Flanges 
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Fig. 10.13. Frequency vs. Displacement for fibre angle rotation in two flanges 


From fig. 10.13, the maximum displacement is observed at natural frequency of 20 Hz. for fibre angle 45° and displacement is 
maximum for fibre angle60° at frequency 5Hz. 
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10.3.2. Effect of fibre angle rotation in web 



Fig. 10.14. Frequency vs. Displacement for fibre angle rotation in web. 
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From fig. 10.14, it is observed that, for fibre angle 45°, the maximum displacement is at natural frequency of 10 Hz. 
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10.3.3. Effect of modular ratio 



Fig. 10.15. Frequency vs. Displacement for various modular ratios. 
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From fig. 10.15. it is observed that, for modulus ratio of 30, the maximum displacement is at natural frequency of 5 Hz. 
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DISCUSSION ON RESULTS 
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DISCUSSION ON RESULTS 

1. It was found that the frequencies of standard beams are higher than that of tapered beams, the frequencies of individual beams 
gradually increases over different sets and also numerical results obtained are almost nearer to the theoretical values obtained 
by using Wittrick.W.H. & Williams,E.W. Algorithm , 

2. The torsional frequencies of first three modes of vibration are increasing significantly for increasing values of taper ratio [] 
and marginally for increasing values of taper ratios a and X. 

3. The torsional frequencies of first three modes of vibration are significantly increasing for increasing values of taper ratios [] 
and X and are decreasing drastically for increasing values of taper ratio a in second mode and third modes of vibration only. 

4. The torsional frequency increases significantly for increasing values of combined influence of negative taper ratios a and [] in 
first mode. In second and third modes of vibrations the torsional natural frequency gradually decreases to certain value and 
then starts increasing. 
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5. An analytical model was developed to study the flexural - torsional vibrations of a laminated composite beam with an I- 
section. The model is capable of predicting accurate natural frequencies as well as vibration mode shapes for various 
configurations including boundary conditions, laminate orientation, ratio of Elastic modulii and height to thickness ratio of the 
composite beams. 
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CONCLUSIONS 
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CONCLUSIONS 


The following conclusions are drawn from the results obtained. 

1. The theoretical results from Finite Element Analysis showed good agreement with experimental results. Therefore, the 
model presented is found to be appropriate and efficient in analyzing free vibration problem of a thin walled laminated 
composite beam. The first four natural frequencies of the rectangular glass epoxy composite beam are found by ANSYS to 
verify with experimental results and the variation is found to be up to 11%. 

2. From the results, it is very clear that changes in laminate sequences yield to different dynamic behaviour of the doubly 
symmetric thin-walled composite I-beams. From this it is concluded that different natural frequencies can be obtained for 
the same geometry, mass and boundary conditions. 

3. The possibility to change fibre orientations, height to thickness ratio and elastic modulus ratio enables to manufacture thin- 
walled composite I-beams to get the desired natural frequencies without increasing its mass and geometry of the beam. 

4. The experimental results and simulation results are found to be good with a variation up to a maximum of 11% .This is due 
to the improper fixing of beam at one end and also may be due to improper use of rubber hammer. 

5. Modal analysis indicated that [0°] s and [0/90] s laminate had higher natural frequencies for vibration and flexural modes 
when compared with the [+/- 45] s laminate, this is due to the higher longitudinal elastic modulus of the laminates. 
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6. The [+/- 45] s laminate had higher torsional modal frequencies due to its higher shear modulus. Harmonic analysis 
confirmed that large increases in normal displacement, rotational displacement, normal stresses and shear stress can occur 
at the resonance frequencies under dynamic loading conditions and these effects must be considered carefully in the design 
mechanical structures that use composite materials. 

7. Harmonic analysis confirmed that large increases in normal displacement, rotational displacement, normal stresses and 
shear stress can occur at the resonance frequencies under dynamic loading conditions and these effects must be considered 
carefully in the design mechanical structures that use composite materials. 

8. It is concluded from the results that [+ 45] laminate is more stable under loading than [0/90] laminate, due to the 
arrangement of the layers. 
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SCOPE OF FUTURE WORK 
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SCOPE OF FUTURE WORK 


• The thickness of web for composite I-beam is considered as constant in the thesis. In future, the mode shapes and its behaviour 
can be studied for a tapered beam when the web thickness tapers uniformly from one end to other end so that weight of the 
beam can be further reduced. 

• the work can be extended for the case of non linear torsional vibration of thin walled composite cantilever I-beam, for both 
symmetric and unsymmetric beams 

• The work can be extended to study the effect of beams resting on continuous elastic foundation or Pasternak foundation or 
generalized foundation on natural frequencies of composite I beams can be studied. 

• The effect of elastic end restraints on the frequency of symmetric and unsymmetric composite I beams can be studied. 
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